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Title: Quantum Multiplexers, Parrondo Games, and Proper Quantization 

A quantum logic gate of particular interest to both electrical engineers and game the- 
orists is the quantum multiplexer. This shared interest is due to the facts that an arbitrary 
quantum logic gate may be expressed, up to arbitrary accuracy, via a circuit consisting 
entirely of variations of the quantum multiplexer, and that certain one player games, the 
history dependent Parrondo games, can be quantized as games via a particular varia- 
tion of the quantum multiplexer. However, to date all such quantizations have lacked a 
certain fundamental game theoretic property. 

The main result in this dissertation is the development of quantizations of history de- 
pendent quantum Parrondo games that satisfy this fundamental game theoretic property. 
Our approach also yields fresh insight as to what should be considered as the proper 
quantum analogue of a classical Markov process and gives the first game theoretic mea- 
sures of multiplexer behavior. 
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Chapter 1 



QUANTUM MECHANICS AND COMPUTATION 



Advances in computation technology over the last two decades have roughly followed 
Moore's Law, which asserts that the number of transistors on a microprocessor doubles 
approximately every two years. Extrapolating this trend, somewhere between the years 
of 2020 and 2030 circuits on a microprocessor will measure on an atomic scale. At 
this scale, quantum mechanical effects will materialize, and virtually every aspect of 
microprocessor design and engineering will be required to account for these effects. 

To this end, quantum information theory studies information processing under a 
quantum mechanical model. One goal of the theory is the development of quantum 
computers with the potential to harness quantum mechanical effects for superior com- 
putational capability. In addition, attention will have to be paid to quantum mechanical 
effects that may obstruct coherent computation. 

The study of possible development of quantum computers falls under the theory 
of quantum computation, an implementation of quantum information theory. Quantum 
computation model quantum information units, called qudits, as elements of a projective 
d-dimensional complex Hilbert space. Physical operations on the qudits are represented 
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by unitary matrices and are viewed as quantum logic circuits. Major results in quan- 
tum computing demonstrate properties of quantum information that are not endemic to 
classical information. Contemporary data implies that in various aspects, quantized in- 
formation offers advantages over classical information. For example, the Deutsch-Jozsa 
quantum algorithm [9] determines whether a function of n binary variables has a spe- 
cific property or not in only one evaluation of the function, compared to the 2"~^ + 1 
evaluations required by the deterministic non-quantum algorithm. Similarly, Grover's 
quantum search algorithm [13] searches a list in time that is quadratic rather than ex- 
ponential in the number of elements in the list, and Shor's period finding quantum al- 
gorithm [32] gives a polynomial time algorithm for factoring integers. The last two 
are well known results in quantum computation and they show that quantum algorithms 
have the potential to out perform classical algorithms for practical problems. 

Quantum game theory offers an exciting and relatively new game theoretic perspec- 
tive on quantum information. Typically, research in the subject looks for different than 
usual behavior of the payoff function of a game when the game is played in a quantum 
mechanical setting. In multi-player games played in a quantum mechanical setting, the 
different than usual behavior of the payoff function studied is typically the occurrence of 
Nash equilibria that are absent in the original game [10, 18, 19]. Because quantum game 
theory has traditionally been heuristic in nature, confusion about and controversy over 
the relevance of "quantum games" to game theory abounds. A resolution to this confu- 
sion and controversy has been recently proposed by Bleiler in [5] via a mathematically 
formal approach to quantum game theory. 

Using Bleiler' s mathematically formal approach to quantum game theory as a step- 
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ping stone, this dissertation promotes the philosophy that quantum game theory should 
be used to gain insights into quantum computation. To this end, the reader is provided 
with a basic introduction to quantum computation and quantum mechanics in the re- 
maining sections of this chapter. In Chapter 2, the Bleiler formalism is reproduced to 
give readers a mathematically formal game-theoretic perspective on quantum games. 
Chapter 3 presents the main results, which are construction and and game theoretic 
analysis of quantum versions of certain one player games, known as history depen- 
dent Parrondo games, and their randomized sequences using the Bleiler formalism as 
a blueprint. These constructions utilize a particular version of a quantum logic cir- 
cuit known as the quantum multiplexer. The connection between quantum game theory 
and quantum computation is made apparent in Chapter 4, where the importance of the 
quantum multiplexer to quantum computation is established via abstract realization of 
an arbitrary quantum logic circuit in terms of circuits composed entirely of quantum 
multiplexers. Chapter 5 may be treated as a stand alone chapter; it proposes the analy- 
sis of quantum circuits acting on exactly two quantum informational units (qubits) via 
quaternionic coordinates. 

1.1 Introduction to Quantum Computation 

Like geometry, quantum mechanics is best viewed axiomatically. For the axioms of and 
basic facts about quantum mechanics, the reader is referred to [26, 27, 6]. These axioms 
and some of the basic facts appear explicitly in the next section during the development 
of one qubit quantum computation. 

A d-ary quantum digit, or qudit for short, is a vector in a complex projective d- 
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dimensional Hilbert space Hd, called the state space of the digit, equipped with the 
orthogonal computational basis 

{|o),|i),...M-i)} 

where |i) = (0, 0, . . . , 1, . . . , 0)^ with a 1 in the (i + l)-st coordinate, for < i < 
{d — 1). To pass from classical to quantum computing, replace a classical rf-ary digit 
(dit) with a qudit as an information unit. The replacement amounts to identifying all 
possible values of the dit with the elements of the computational basis of the state space 
of the corresponding qudit. This identification enlarges the set of operations on the 
dit to include quantum operations which, by the axioms of quantum mechanics, are 
represented by unitary operators on the state space. One then typically explores whether 
this enlargement results in any computational advantages or enhancements. 

To be more specific, unitary operators can be used to create complex projective linear 
combinations of the basis qudits. In other words, a qudit \a) in Hd can be expressed as 
a complex projective linear combination of the basis qudits 

\a) = XiEC 

1=0 

where \a) = X\a) for any non-zero complex number A. Physicists call this complex 
number A a phase. Up to phase, the state \a) can be normalized; that is, \a) can be 
expresses with 

d-l 
1=0 
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The measurement axioms of quantum mechanics say that the real number is the 
probabiUty that the state vector | a) will be observed in i-th basis state upon measurement 
with respect to that basis. Typical considerations in quantum computing are whether 
evolutions of the state space offer computational enhancements. 

When considering several qudits at once, the axioms of quantum mechanics tell us 
to consider their joint state space. When the state spaces of n qudits of different d- 
valued dimensions are combined, they do so via their tensor product as per the axioms 
of quantum mechanics and the result is a n qudit hybrid state space 

n^Hd^^Hd,®---® Hd^ 

where 7^^. is the state space of the -valued qudit. The computational basis for Ti, con- 
sists of all possible tensor products of the computational basis vectors of the component 
state spaces Hdi- If di — d for each i, the resulting state space TYf" is that of n d- valued 
qudits. 

Once a basis for the state space has been chosen, a unitary operator on it is repre- 
sented by a unitary matrix. For the hybrid state space H, an evolution matrix will be 

of size {did2 ■ ■ ■ c^at) x {did2 ■ ■ ■ d^), while the evolution matrix for 7Y|" will have size 
X d". 

Consider a two dimensional state space H2- This is the state space of a quantum 
system which gives two possible outcomes upon measurement. An example of such 
a system would be one that describes the spin states of an electron. Topologically, 
H2 — CP^. The two possible states of the system form the computational basis for the 
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state space. These orthogonal basis state are viewed as the two possible values a bit of 
information can take on. Call the elements of 7^2 qubits, short for binary quantum digit. 
The resulting 2-valued quantum computing has traditionally been the most active area 
of research. The basics of 2-valued quantum computing are reviewed in the following 
sections. Higher valued quantum computing has seen much research activity recently as 
well. The reader is referred to chapter 2 for a discussion of certain aspects of valued 
quantum computing and relevant references. 

1.1.1 One Qubit Quantum Computing 

Let |6o) and be an orthogonal basis for ^,2- Then the states of the qubit are projective 
linear combinations of these basis elements over C: 



1 2 1 2 

withao)<^i ^ C satisfying, without loss of generality, |q;o| +|cti| = 1. These projective 
complex linear combinations are also called superpositions of the states |0) and 1 1) . The 
computational basis is the set 



comp 



6 



Chapter 1. Quantum Mechanics and Computation 



which gives the convention of labeling the basis with Boolean names, with 



\t>o) = |0) 



and 1 61) = |1) 



V V 



But note that these are only names. For example, in the spin state model for an 
electron, one might imagine that |0) is being represented by an up-spin while |1) by a 
down-spin. The key is that there is an abstraction between the technology (spin state 
or other quantum phenomena) and the logical meaning. This same detachment is true 
in classical computers where we traditionally call a high positive voltage "1" and a low 
ground potential "0". 

Let iV'i) = ao |0) + ai\l) and 



U 



—r] 



\ 



7] 



(1.1) 



be a special unitary operator which, by axioms of quantum mechanics, corresponds to a 
physical operation. Further, suppose that 77 is a complex root of unity other than ±1, the 
use of which will be justified shortly. The matrix U acts on |^i) as follows. 



-f] 
77 



\ / \ / - \ 

-riai 



I 



ttO 



770:0 



-77q;i |0) +77q;o |1) 



(1.2) 



Up to multiplication by unitary phase, the operator U interchanges the coefficients of 
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Figure 1.1: Inverter or the NOT quantum logic gate U. The wires carry quantum information, 
namely qubits. 





|0, 

Figure 1.2: Standard notation for the NOT gate, 
the basis states of CP^. In particular, U sends the state |0) to the state rj |1) 

/ 

U\0) = 

\ 

and the state 1 1) to the state —f] |0) 

(I —n 



-f] 

7] 



r/|l) 



U\l) = 




-7] 


u 




1 


-7] 






' 1 


V 


I 







This action of U is interpreted as that of a quantum logic gate that inverts, up to 
unitary phase, the logical values |0) and that is, the gate ^7 is a quantum version 
of the NOT gate in classical logic. This point of view allows one to view quantum 
mechanics as a theory of quantum computation. Standard notation for the NOT gate is 
given in Figure 1.2. 

In the quantum theory of games, one frequently views a qubit as a "quantum coin" 
and hence the gate U can be interpreted as the quantum mechanical analog of flipping 
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over a coin, while the 2x2 identity matrix is the analog of leaving the coin un-flipped. 
In certain quantum games, such as the ones found in [18, 1], the flipping and un-flipping 
actions of players on the so-called maximally entangled state of two qubits 



^(|00) + |11)) 



are considered. For the purpose of analysis of the quantum game, these actions are 
required to produce an orthogonal basis of the joint state space, and this happens only 
when T] is an appropriate root of unity other than ±1. 

1. 1 .2 The One Qubit Hadamard Quantum Logic Gate 

Quantum computing literature gives many interesting examples of one qubit gates. The 
focus here will be on the one qubit Hadamard gate described by the special unitary 
matrix 



V2 



1 1 



Application of the gate H to either basis states |0) and 1 1) creates an equal superposition 
of the basis state, that is, a superposition that will appear in each basis state with equal 
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^(|o>-|i» 




0)+aj 1 



Figure 1.3: The Hadamard gate H that puts a basis state into an equal superposition of the basis 
states, and an arbitrary state \a) into the superposition ao |0) + ai 



probability upon measurement with respect to the basis. 



H\0) 



i 



^ 1 \iA 



1 -1 



V5 



+ 



VI 



and 



H\l) 



I 
i 



^ 1 \ (,^ 

1 -1 



1 



\-v 



Vl 



VI 
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1.1.3 Measurement 

In Equation (1.2), if both ao, ai ^ 0, then the how does one interpret the complex pro- 
jective linear combination — f]ai \ 0) +r]ao 1 1 ) of the basis states in the context of comput- 
ing? The answer comes from quantum mechanics' axiom of measurement which allows 
a probabilistic interpretation of such complex projective linear combinations as follows. 
Upon measurement with respect to the orthogonal basis {|0) , |1)}, the combination is 
observed to be in the basis state |0) with probability \aif and in the basis state 1 1) with 

2 I 1 2 I 1 2 

probability |a;i| (remember that ?7 is a unit complex number so |?7| = |r/| = 1). Two 
important measurement operators are 



Mn 



^ 0^ 



M, 



The measurement operator Mq projects a complex projective linear combination onto 
the basis state |0) while Mi projects onto the basis state For example, let 

1^) = ao|0) +«! |1) 



Then the probability of measuring the complex projective linear combination 
basis state |0) is 



in the 



Pirn 



a b 



a 



1 
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Note that measurement operators are not quantum logic gates as they are non-unitary, 
but rather are projections onto the basis states. 

1.2 Quantum Computing with Multiple Qubits 

Quantum computing can be extended to multiple qubits via the creation of composite 
state spaces from the state spaces of many individual qubits. 

For example, consider two qubits = a |0) + 6 |1) and \ip2) = c |0) + o? both 
written with respect to the computational basis. Then the joint state of the total system 
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is given by: 



iV'i) ^ 1^2) = = ac |0) ® |0) + ad |0) ® |1) + 6c |1) ® |0) + 6d |1) ® |1) 



= ac 



+bd 



' ^ 



' ^ 



+ be 



ac 








+ ad 



1 




6c 





1 

voy 



6rf 






viy 



ac |00) + ad |01) + be |10) + bd\ll) 



1.2.1 Two Qubit Quantum Gates 



An easy way to obtain two qubit quantum gates is by producing the tensor product of 
two one qubit gates. That is, if Ui and U2 are one qubit gates, then 



U^Ui®U2 



is a two qubit gate. Two qubit gates such as U above that are tensor products of one qubit 
gates act locally on each qubit due to the bi-linearity of the tensor product. Nonetheless, 

13 



Chapter 1 . Quantum Mechanics and Computation 





Figure 1.4: The two qubit Hadamard gate is just the tensor product of the one qubit Hadamard 
gates acting on each qubit. In general, multiqubits gates can be created via the tensor product 
of one qubit gates. However, it is not always true that a multiqubit gate is equal to the tensor 
product of one qubit gates. Consider for example the CNOT gate of Figure 1.5. 

such gates are crucial to quantum computing. For example, the two qubit Hadamard 
gate defined as 



V2 
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1 -1 
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V 


1 






V 


-1 ^ 


' I 




V 


1 -1 


-1 1 / 



is essential for the creation of a particular equal superposition of two qubits which plays 
a crucial role in the development of quantum algorithms that out-perform classical al- 
gorithms [13, 32]. 



14 



Chapter 1. Quantum Mechanics and Computation 



1.2.2 Controlled NOT (CNOT) gate 

Perhaps the most important two qubit gate is the controlled NOT (CNOT) gate. Its 
importance lies in its property of forming, together with one qubit gates, sets of universal 
quantum logic gates. Informally, a set of quantum logic gates is universal if any quantum 
logic gate may be approximated by the gates in the set to arbitrary accuracy. For a 
detailed discussion of universality, the reader is refer 

The CNOT gate acts as a NOT gate on the second qubit (target qubit) if the first 
qubit (control qubit) is in the computational basis state |1). So when passing through 
the gate the states 1 00) and 1 01) are unaltered, while the state 1 10)is sent to 1 1 1) and vice 
versa. In the joint computational basis, the CNOT gate is 



CNOT 



/ 1 ^ 

10 

1 

\ 1 y 



Note that the CNOT gate is not the tensor product of any pair of one qubit gates. Indeed, 
there are plenty of other two and multiqubit gates that are not tensor products of one 
qubit gates. This property of quantum logic gates is one more reason that quantum logic 
circuit synthesis is a much studied subject. 
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Figure 1.5: The controlled NOT (CNOT) gate. The vectors |00) and |01) are unaltered, while 
the vector |10) is sent to |11) and vice versa. 

1.2.3 Entanglement 

Entanglement is a uniquely quantum phenomenon. Entanglement is a property of a 
multi-qubit system and can be thought of as a resource. To explain entanglement, let us 
examine a so-called EPR pair of qubits named after Einstein, Podolsky, and Rosen. The 
CNOT gate will be used in this example. 

We begin with two qubits j^'i) = |0) and IV'2) = Apply the Hadamard gate to 
1-01) to get 

W.) = H\^l^,) = ^^ + ^\l) 
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H 











^(loo)^")) 




Figure 1.6: The Hadamard gate H that puts a basis state into an equal superposition of the basis 
states, and an arbitrary state |a) into the superposition ao |0) + ai |1). 



The joint state-space vector is the tensor product 



IV'i) ® \i^2) = WM = ^ |oo) + (0) |oi) + |io) + (0) 111) 



Now apply the CNOT gate to this joint state of the two qubits. This gives 
/ 



1 ^ 

10 

1 

\ 1 y 



/ J- \ 





V2 



( J- \ 




V 71 / 



I 

72 



(l\ 






75 








The final joint state above has the property that it cannot be built up from the tensor 
product of states in the component spaces of each qubit. That is. 



-J=|00) + -J=|11)^ |0l)®|02). 
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To illustrate why entanglement is so strange, let's consider performing a measure- 
ment just prior to applying the CNOT gate. The two measurement operators (for obtain- 
ing a |00)ora|ll)) are: 



Moo = 



/ 1 ^ 



10 

y y 



and Mil = 





10 



y 1 y 



Just prior to the CNOT the system is in the state 



therefore 



-^|00) + 0|01) + -^|10)+0|11), 



p(0) 



— — 



1 ^ 



10 

\ y 







V2 
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Hence the result of measuring will clearly be |0). After the measurement, we have 







V2 



and we see that measurement had no effect on the first qubit and it remains in a super- 
position of |0) and |1). Now consider the same measurement but just after the CNOT 
gate is applied, with the joint state iV's) = |00) + 



p(0) 



— n n — 



1 

10 





1 ^ \ 















1 







~ 2 




i 





V2 



Hence, after the CNOT gate is applied we have only a 50% chance of obtaining |0). Of 
particular interest to our discussion, however, is what happens to the state vector of the 
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system after measurement. 













i|00) 



This is the remarkable thing about entanglement. By measuring one qubit we can affect 
the probability of the state observations of the other qubits in a system! The state of the 
other qubit |^^) = -^^ |0) + |1) is changed to |0) after the measurement. 
Quoting Oskin [27] regarding entanglement: 

"How to think about this process (entanglement) in an abstract way is an 
open challenge in quantum computing. The difficulty is the lack of any 
classical analog. One useful, but imprecise way to think about entangle- 
ment, superposition and measurement is that superposition "is" quantum 
information. Entanglement links that information across quantum bits, but 
does not create any more of it. Measurement "destroys" quantum informa- 
tion turning it into classical. Thus think of an EPR pair as having as much 
"superposition" as an un-entangled set of qubits, one in a superposition be- 
tween zero and one, and another in a pure state. The superposition in the 
EPR pair is simply linked across qubits instead of being isolated in one." 
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A FORMAL APPROACH TO QUANTUM GAMES 



One way to view a game is as a function. We view here quantum games as extensions 
of such functions. For a detailed and formal introduction to game theory the reader 
is referred to [ ] and [_-r]. The following discussion on quantum games that follows is 
motivated by a mathematical formalism for "quantum mixtures" developed by S. Bleiler 
in [. ] and reproduced in section 2.1 below. 

Recall that a key goal in the study of multi-player, non-cooperative games is the 
identification of potential Nash equilibria. Informally, a Nash equilibrium occurs when 
each player chooses to play a strategy that is a best reply to the choice of strategies of 
all the other players. In other words, unilateral deviation from the choice of strategy at a 
Nash equilibrium by any player is detrimental to that player's payoff in the game. How- 
ever, in finite classical games, Nash equilibria may not exist. In such situations, classical 
game theory calls upon the players to randomize between their strategic choices, also 
known as mixing strategies. For finite games, Nash proved [25] that this gives rise to 
Nash equilibria in the "mixed game" that simply do not exist in the original game. For- 
mally, the mixed game is the result of an extension of the payoff function of the original 
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game to a larger set of strategies for each player. 

The Bleiler formalism for quantum mixtures views quantum game theory in this 
light. That is, this formalism views quantum game theory as an exercise in the extension 
of the payoff function of a game with the goal of finding Nash equilibria with higher 
payoffs that were un-attainable in the original game or its "classical extensions". The 
extensions dealt with in quantum game theory are referred to as a quantization pro- 
tocols. This mathematically formal perspective provides a game theoretic context in 
which many issues in quantum game theory can be discussed and potentially resolved. 
For example, critics of quantum game theory wonder whether instances of Nash equi- 
Ubira with higher payoffs in certain quantum games are just Nash equilibria of some 
other classical game theoretic construction realized quantum mechanically. This point 
of view implies that quantum game theory is essentially a study in expensive ways to 
generate classical game theoretic results and offers nothing "new" to game theory. 

Such criticism is addressed in the Bleiler formalism which points out that any quan- 
tum game that contains the original or the classical game as an embedded subgame has 
the potential to offer something new to the game's analysis. When a quantum game has 
this property, it is referred to in the Bleiler formalism as a proper quantization of the 
original game. When a quantum game carries an embedded copy of the mixed version 
of the original game, the formalism refers to it as a complete quantization of the orig- 
inal game. Much of the current work in quantum game theory can be characterized as 
calUng upon the players to use the higher orders of randomization given by quantum 
superpositions and randomized quantum superpositions. Call these quantum strategies 
and mixed quantum strategies, respectively. If the quantization of the game is proper 
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or complete, then any new Nash equilibria with higher payoffs that result from the use 
of quantum or mixed quantum strategies can be meaningfully compared with the Nash 
equilibria of the original game. 

A detailed review of the Bleiler formalism follows. 

2.1 The Bleiler Formalism for Quantum Mixtures 

Definition 2.1. Given a set {1,2, • • • ,n} of players, for each player a set Si {i = 
1, • • • , n) of so-called pure strategies, and a set ilj (i = 1, • • • , n) of possible outcomes, 
a game G is a vector- valued function whose domain is the Cartesian product of the Si's 
and whose range is the Cartesian product of the ili's. In symbols 

n n 

G.HSi^Hni 

i=l 1=1 

The function G is sometimes referred to as the payoff function. 

Here a play of the game is a choice by each player of a particular strategy Si the 
collection of which forms a strategy profile (si, • • • , s„) whose corresponding outcome 
profile is G{si, ■ ■ ■ ,Sn) — ((^i, • • • , ojn), where the UiS represent each player's individ- 
ual outcome. Note that by assigning a real valued utility to each player which quantifies 
that player's preferences over the various outcomes, we can without loss of generality, 
assume that the Qj's are all copies of R, the field of real numbers. 

In game theory, players' concern is the identification of a strategy that guarantees a 
maximal utility. For a fixed (n — 1) -tuple of opponents' strategies, rational players seek 
a best reply, that is a strategy s* that delivers a utiHty at least as great, if not greater, 
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than any other strategy s. When every player can identify such a strategy, the resulting 
strategy profile is called a Nash equilibrium. Formally, 

Definition 2.2. Let s_j be a strategy profile of all players except player i. A Nash 
equilibrium (NE) for the game G is a strategy profile (s*, S-i) such that 

where for all i, Si,s* G Si and s* ^ Si. 

Other ways of expressing this concept include the observation that no player can 
increase his or her payoffs by unilaterally deviating from his or her equilibrium strategy, 
or that at equilibrium all of a player's opponents are indifferent to that player's strategic 
choice. As an example, consider the Prisoner's Dilemma, a two player game where each 
player has exactly two strategies (a so-called 2 x 2 or bimatrix game) and whose payoff 
function is indicated in Table 2.1. The rows of Table 2.1 contain the strategies of player 
1 while the columns contain the strategies of player 2. 

Note that for player 1 the pure strategy S2 always delivers a higher outcome than 
the strategy si (say $2 strongly dominates si) and for player 2 the strategy t2 strongly 
dominates ti. Hence the pair (52,^2) is a (unique) Nash Equilibrium. 

However, games need not have equilibria amongst the pure strategy profiles as ex- 
emplified by the 2x2 game of Simplified Poker whose payoff function is given in Table 
2.2. 

As remarked above, the game theoretic formalism now calls upon the theorist to 
extend the game G by enlarging the domain and extending the payoff function. Of 
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h 


t2 






Sl 


(3,3) 


(0,5) 








(5,0) 


(1,1) 




Table 2.1: 


Prisoner's Dilemma 




h 


t2 




(5/4, 


-5/4) 


(0,0) 


S2 


(0,0) 


(5/2, - 


-5/2) 



Table 2.2: Simplified Poker. 



course, the question of if and how a given function extends is a time honored problem 
in mathematics and the careful application of the mathematics of extension is what will 
drive the formalism for quantization. Returning to classical game theory, a standard 
extension at this point is to consider for each player the set of mixed strategies. 

Definition 2.3. A mixed strategy for player i is an element of the set of probability 
distributions over the set of pure strategies Si. 

For a given set X, denote the probability distributions over X by A(X) and note 
that when X is finite, with k elements say, the set A(X) is just the k — 1 dimensional 
simplex A*^*^~^^ over X, i.e., the set of real convex linear combinations of elements of 
X. Of course, we can embed X into A.(X) by considering the element x as mapped to 
the probability distribution which assigns 1 to a; and to everything else. For a given 
game G, denote this embedding of Si into A(5'i) by e^. 

Let p = (pi, . . . ,p„) be a mixed strategy profile. Then p induces the product dis- 
tribution over the product J| Si. Taking the push out by G of the product distribution 
(i.e., given a probabiHty distribution over strategy profiles, replace the profiles with their 
images under G) then gives a probability distribution over the image of G, ImG. Fol- 
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Product 



Push-out 



^ A(lmG) 




E 



G 



Figure 2.1: Extension of the game G to G' 



"mtx 



lowing this by the expectation operator E, we obtain the expected outcome of p. Now 
our game G can be extended to a new, larger game G^^^ . 

Definition 2.4. Assigning the expected outcome to each mixed strategy profile we obtain 
the extended game 



Note G*™*^ is a true extension of G as G*™*^ o Ilej = G\ that is, the diagram in Figure 
2.1 is commutative. 

As remarked above, Nash's famous theorem [25] says that if the Si are all finite, then 
there always exists an equilibrium in G""*"^. Unfortunately, this equilibrium is called a 
mixed strategy equilibrium for G, when it is not an equilibrium of G at all, the abusive 
terminology confusing G with its image, ImG. 

2.1.1 Quantization 

The Bleiler formalism asserts that some of the controversies surrounding quantum game 
theory may be resolved if one focuses on the quantization of the payoffs of the original 
game G, and expresses the quantized version of G as a (proper) extension of the original 
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payout function in the set-theoretic sense, just as in the classical case. 

Classically, probability distributions over the outcomes of a game G were con- 
structed. Now the goal is to pass to a more general notion of randomization, that of 
quantum superposition. Begin then with a Hilbert space Ti, that is a complex vector 
space equipped with an inner product. For the purpose here assume that H is finite di- 
mensional, and that there exists a finite set X which is in one-to-one correspondence 
with an orthogonal basis B of H. When the context is clear as to the basis to which 
the set X is identified, denote the set of quantum superpositions for X as QS{X). Of 
course, it is also possible to define quantum superpositions for infinite sets, but for the 
purpose here, one need not be so general. What follows can be easily generaUzed to the 
infinite case. 

As mentioned above, the underlying space of complex linear combinations is a 
Hilbert space; therefore, we can assign a length to each quantum superposition and, up 
to phase, always represent a given quantum superposition by another that has length 1. 

For each quantum superposition of X we can obtain a probability distribution over 
X by assigning to each component the ratio of the square of the length of its coefficient 
to the square of the length of the combination. This assignment is in fact functional, and 
is abusively referred to as measurement. Formally: 

Definition 2.5. Quantum measurement with respect to X is the function 

q~ : QS{X) — > A{X) 
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given by 



ax + (3y 




) 



Note that geometrically, quantum measurement is defined by projecting a normal- 
ized quantum superposition onto the various elements of the normalized basis B. Denote 
quantum measurement by q'^^'^'^ if the set X is clear from the context. 

Now given a finite n-player game G, suppose we have a collection Qi, . . . , Q„ of 
non-empty sets and a protocol, that is, a function 6 : H 2* ^ QS'(ImG). Quantum 
measurement q^^Q then gives a probability distribution over ImG. Just as in the mixed 
strategy case we can then form a new game by applying the expectation operator E. 

Definition 2.6. Assigning the expected outcome to each probability distribution over 
ImG that results from quantum measurement, we obtain the quantized game 



Call the game G® thus defined to be the quantization of G by the protocol Q. Call 
the Qj's sets of pure quantum strategies for G®. Moreover, if there exist embeddings 
e'j : Si — > Qi such that G® o H — ^ proper quantization of G. If there exist 

embeddings e^' : A{Si) Qi such that G® o [] e'l = G™^, call G® a complete quan- 
tization of G. These definitions are summed up in the commutative diagram of Figure 
2.2. Note that for proper quantizations, the original game is obtained by restricting the 
quantization to the image of Y[ e';. For general extensions, the Game Theory literature 
refers to this as "recovering" the game G. 

It follows from the definitions of G™^ and G® that a complete quantization is proper. 
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QS(lmG) 




Figure 2.2: Extension of the game G to G®. 



Furthermore, note that finding a mathematically proper quantization of a game G is 
now just a typical problem of extending a function. It is also worth noting here that 
nothing prohibits us from having a quantized game G® play the role of G in the classical 
situation and by considering the probability distributions over the Qi, creating a yet 
larger game C"®, the mixed quantization of G with respect to the protocol 0. For a 
proper quantization of G, G""® is an even larger extension of G. The game G""® is 
described in the commutative diagram of Figure 2.3. 

In many cases, the Qi of the quantization protocols are expressed as quantum opera- 
tions. These operations require a state to "operate" on. In this situation the definition of 
protocol additionally requires the definition of an "initial state" together with the fam- 
ily of quantum operations which act upon this state, along with a specific definition of 
how these quantum operations are to act. As exemplified in the next chapter, different 
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Figure 2.3: Extension of the game G® to G' 



{Q,I) ■ ^^(imG) ^^"^ A (Im G ) 




Figure 2.4: Proper quantization of a one player game with strategy space S via the protocol 6 
and quantum sttategy space Q. 

choices for the initial state can give rise to very different protocols sharing a common 
selection and action of quantum operations. When a protocol depends on a specific 
initial state /, the protocol is then denoted by 6/. 

In subsequent sections, a version of this formalism adapted to one player games will 
be utilized. The underlying quantization paradigm being the replacement of probability 
distributions by the more general notion of quantum superposition followed by mea- 
surement. The functional diagram for proper quantization that will be utilized is given 
in Figure 2.4 where the commutativity of the diagram requires that E o (gj^^^) o 6 o e = 

o e = G. Incorporating the discussion above, when games Gs and protocols 0/ de- 
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pend on a given initial states s and /, respectively, the initial states s and / are regarded 
as part of the single player's strategic choice. In these cases, the embedding e of 5" into 
Q additionally requires the mapping of the initial state s of Gg to the initial state / of 
the protocol ©7. The resulting quantum game is denoted by Gf'. 
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PROPERLY QUANTIZING HISTORY DEPENDENT PARRONDO GAMES 



A major insight about quantized games that results from the Bleiler formalism discussed 
in Chapter 2 is that for the quantization of a game to be game-theoretically significant, 
it must be proper. Previous work on the quantization of the history dependent Parrondo 
game by Flitney, Ng, and Abbott (FNA) [11] produced quantizations that are not proper. 
In this chapter, after recalling the basic facts regarding Parrondo games and the FNA 
quantization protocols, proper quantizations for the history dependent Parrondo game 
and their randomized sequences are constructed. 

3.1 Parrondo Games 

Parrondo et. al first formulated such games in The subject of Parrondo games has 
seen much research activity since then. Parrondo games typically involve the flipping of 
biased coins and yield only expected payoffs. A Parrondo game whose expected payoff 
is positive is said to be winning. If the expected payoff is negative, the game is said to 
be losing, and if the expected payoff is 0, the game is said to he fair. 

Parrondo games are of interest because sequences of such games occasionally ex- 
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hibit the Parrondo effect; that is, when two or more losing games are appropriately 
sequenced, the resulting combined game is winning. Frequently, this sequence is ran- 
domized which means that the game played at each stage of the sequence is chosen 
at random with respect to a particular probability distribution over the games being se- 
quenced. A comprehensive survey of Parrondo games and the Parrondo effect by Abbott 
and Harmer can be found in [ ]. 

Earlier work on the quantization of Parrondo games can be found in ' ] where 
Meyer offers an analysis of a quantization of a particular type of Parrondo game, and in 
[11] where Abbott, Flitney, and Ng (AFN) propose quantizations of a different type of 
Parrondo game. The authors of both papers quantize their original game via their own 
particular quantization protocols, and further, model the game sequences as iterations 
of their protocols. In each of these protocols, quantum actions are performed on a 
collection of initial states of a quantum system. At the end, a measurement of certain 
specific states is made and, from the resulting probability distributions, an expected 
payoff computed. 

3.1.1 Capital Dependent Parrondo Games 

In [29], Parrondo et al describe two types of coin flipping games which have the property 
that if individually repeated, the games result in a decreasing expected payoff to the 
player, yet when the two games are played in a deterministic or probabilistic sequence 
repeatedly, the expected payoff to the player increases over time. 

Suppose that X(t) = 0, 1, 2, . . . is the capital available to the player. If the player 
wins a game, then the capital increases by one, and if the player loses, then the capital 
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Table 3.1: Game B 





Prob. of gain 


Prob. of loss 


X{t) = mod 3 


Pi 


1 -pi 


X{t) = 1 or 2 mod 3 


P2 


1 -P2 



decreases by one. The simplest type of this game, referred to in the literature as game A, 
is determined by a biased coin with probability of gain p. That is, the capital increases 
by one with probability p and decreases by one with probability 1 — p. Another game, 
called game B, is defined by two biased coins. The choice of which coin is to be played 
in an instance of the game B is determined by the congruence modulo 3 of the capital, 
X{t), available to the player in that instance. Hence, game B is defined by the rules 
given in table 3.1. 

Parrondo et al set p = ^ — e, j9i = ^ — e, j92 = | — e, for e > as an example 
of games A and B which are losing if played individually or in a fixed sequence, but 
which, when combined in a randomized sequence with the uniform distribution over the 
two games, is winning. Both games A and B are losing, winning, or fair as e > 0, e < 
and e = 0, respectively. Parrondo et al consider in detail the case when both games A 
and B are fair. The game B is analyzed as a Markov process Y{t) = X{t) (mod 3), 
that is, Y{t) is equal to the remainder upon dividing the capital X(t) by 3. A transition 
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matrix for game B is thus given by 



T 



^ 1 5\ 

"44 



^ i 
10 " 4 



9 3 n 

V To 4 / 



(3.1) 



The stationary state for this Markov process can be computed from the matrix equation 



^ ^ 3 \ / \ / \ 

U 4 4 TTo TTq 



^ i 
10 " 4 



9 3 



TTl 



TTl 



(3.2) 



Vto 4 uyyTTsy j 



where vTj is the probability of the capital X{t) taking on a value congruent to i (mod 3), 
i = 0, 1, 2. The matrix Equation (3.2) gives rise to the following system of equations 



-vri + -7r2 = vro 

9 3 

TTn H TTi = TTo 

10 4 



(3.3) 
(3.4) 
(3.5) 
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which has the following solution. 



TTo = TTo, TTi = -TTq, 772 = ^TTq. 



Since the game is assumed to be fair, pi7ro+p27ri+p27r2 = |, and one computes ttq 



_5_ 
13' 



Now if the fair games A and B are played in a randomized sequence, the resulting 
capital can be increasing. To see this, let q be the probability with which the game A is 
played. Then game B is played with probability (1 — q). Again, analyze the Markov 
sequence Y{t) = X{t) (mod 3), but this time the transition matrix is 



/ 



r 







\ 



(3.6) 



To sequence these games via the uniform distribution, set q = \ and get 



r 



4 8 4 8 ^ 



4 ~ 20 



1 _|_ _9_ 1 _|_ 3 
\ 4 ~'~ 20 4 ~'~ 8 



10 ^ 8 



7 5 n 

V 10 8 U / 



(3.7) 
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Before last 
t-2 



Last 

t-1 



Coin 



Prob. of gain 
att 



Prob. of loss 
at t 



gain 
gain 
loss 
loss 



gain 
loss 
gain 
loss 




Pi 

P2 
P3 
Pi 



1 -pi 
1-P2 
1 -P3 
1-P4 



Table 3.2 



History dependent game B'. 



Computing the stationary state (ttq, n[, Tx'2f for the case in which each game A and 
B is fair, gives ttq = vr'^ = and 'n'2 = up to a normalization constant. Note 
that ttq = 1^ is larger than ^, and thus the capital increases. 

3.1.2 A History Dependent Parrondo Game 

The history dependent Parrondo game, introduced in \^'^] by Parrondo et al, is again 
a biased coin flipping game, where now the choice of the biased coin depends on the 
history of the game thus far, as opposed to the modular value of the capital. A history 
dependent Parrondo game B' with a two stage history is reproduced in Table 3.2. 

As above, let X{t) be the capital available to the player at time t. At stage t, this 
capital goes up or down by one unit, the probability of gain determined by the biased 
coin used at that stage. Obtain a Markov process by setting 



This allows one to analyze the long term behavior of the capital in game B' via the 




(3.8) 
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stationary state of the process Y{t). The transition matrix for this process is 



X 



Pi 



P2 



P3 

1-P3 

PA 





\ 1-P2 I -Pa J 
The stationary state can be computed from the following equations 



(3.9) 



and is given by 



(1 - pi)7ri + (1 - Pz)'Kz = TTa 

P27r2 + PaT^A = TTs 
(1 - P2)'K2 + (1 - Pa)t^A = 







^ PsPa ^ 




1 


P4(l -Pl) 




~ TV 


P4(l -Pl) 






(i-pi)(i-P2) y 



(3.10) 



after setting the free variable va — {1 — pi){l — P2) and normalization constant 



N 



. Y.{^0? = V (^^3P4)2 + 2 [(1 - p,)pAf + [(1 - Pl) (1 - P2t 

\ i=l 
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which simpUfies to 

AT = (1 - pi) (2p4 + 1 - P2) + P3P4- 

Consequently, the probabiUty of gain in a generic run of the game B' is 

B' - Jl - P4 (P3 + 1 - Pi) , , . 

P.^^ - Nf-. ""^P^ - (1 - pO (2p4 + 1 - P2) + P3P4 ^ ^ 

where tt^ is the probability that a certain history j, represented in binary format, will 
occur, while pj is the probability of gain upon the flip of the last coin corresponding to 
history j. The expression for p^^^ simplifies to 

Pfa'in = 1/(2 + a;/!/) (3.12) 

with 

y = P4(P3 + l-Pi) > (3.13) 
for any choice of the probabiUties pi, . . . P4, and 

X = (1 - pi)(l - P2) - P3P4. (3.14) 

Therefore, game B' obeys the following rule: if x < 0, B' is winning, that is, has 
positive expected payoff; if a; = 0, B' is fair; and if a; > 0, B' is losing, that is, has 
negative expected payoff. 

Before proceeding further, it is useful to view the preceding ideas in a more formal 
game theoretic context. For this, consider the Parrondo games as one player games in 
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normal form, that is, as a function, where the one player's strategic choices in part cor- 
respond to the biases of the coins. For a history dependent Parrondo game with two 
historical stages, Parrondo et al refer to these choices as a "choice of rules." However, 
the mere choice of biases for the coins is not enough to determine a unique normal form 
for these history dependent Parrondo games. In particular, an initial probability distri- 
bution over the allowable histories is also required. Although any specific distribution 
suffices to uniquely determine such a normal form, as the structure of the game is given 
by a Markov process, there is a natural choice for this initial distribution. Though this 
issue is not discussed by Parrondo et al, these authors immediately focus on this natu- 
ral choice, namely, the distribution corresponding to the stationary state of the Markov 
process representing the game. 

Now, the normal form of these history dependent Parrondo games maps the tuple 
(P, s) into the element 

(tTiPi, 7ri(l - pi), 7I2P2, Vr2(l - P2), TTsPs, 7r3(l - Pa), TTiP^, TTi{l - Pi)) 

of the probability payoff space [0, 1] where s = (tti, 7r2, tts, 114) E A(hist(j') is the sta- 
tionary state of the Markov process with transition matrix defined by P = {pi,P2, Pa, P4) 
e [0, 1]^"^, as in Equation (3.9). Formally, 

a : [0, 1]^^ X A(histG') ^ [0, 1]^^ (3.15) 

Gs : (P, s) ^^ (7ripi,7ri(l - pi), 712^2, vr2(l - P2), TTspa, 713(1 - ps), 7r4p4, 7r4(l -^4)) 

(3.16) 
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The outcomes winning, breaking even, or losing to the player occur when p^^^^ > |, 
PgL = h andpfa'in < |, respectively. 

Note that in this more formal game theoretic context for history dependent Parrondo 
games, the dependence of these games on the initial probability distribution s is made 
clear. This initial probability distribution plays the role of the initial state s for the 
classical game Gg appearing in the proper quantization discussion at the end of chapter 
2. 



3.1.3 Randomized Combinations of History Dependent Parrondo Games 

Consider now the two stage history dependent game obtained by randomly sequencing 
the games B' and B" where each of B' and B" are history dependent Parrondo games 
with two stage histories. This can be formally considered as a real convex linear com- 
bination of the games B' and B", where the coefficients on B' and B" are given by r, 
the probability that the game B' is played at a given stage, and (1 — r), the probability 
that the game B" is played at a given stage. This is because the transition matrix of the 
Markov process associated to the randomized sequence is obtained from the transition 
matrices T' and T" for the games B' and B", respectively, by taking the real convex 
combination rT' + (1 - r)T". Explicitly, let 



T' 



ai 
1 — ai 







Oi2 



0(3 
1-^3 

0:4 



\ l-a2 







(3.17) 



1 — J 
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and 



rpll 



I3i 
1-A 
(52 



I3z 
l-/?3 

/?4 



(3.18) 



\ 1-132 1 - y 

with cXjiPj E [0,1] representing the probability of gain for the j coin in games B' and 
B" respectively. Then the transition matrix rT' + (1 — r)T" of the Markov process for 
the randomized sequence of B' and B" consists of entries tj = raj + (1 — r)(Pj) and 



[1 — aj) + (1 — r)(l — Pj) in the appropriate locations. Call this randomized 



sequence of games B' and B" the history dependent game B'B" with probability of 
gain tj . The stable state, computed in exactly the same fashion as the stable state for the 
game B' in section 3.1.2 above, has form 







/ 
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r2 


1 




U{l-ti) 






~ R 








V n ) 




\ 




/ 



(3.19) 



with R = X]t=i '0 ^ normalization constant. Using the stable state, the probability of 



gain in the game B'B" is computed to be 

-± 



B'B" 
gain 



{l-t,){2U + l-t2)+hU' 



(3.20) 
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Just as in case of the game B', the expression for Pgajn reduces to 

p!2" = 1/(2 + x'/y') (3.21) 

with 

y' = U{ts + 1 - ti) > (3.22) 
for any choice of the probabilities ti, . . .t^, and 

x' = {l-h){l-t2)-tsU. (3.23) 

The game B'B" therefore behaves entirely like the game B', following the rule: if 
x' < 0, B'B" is winning, that is, has positive expected payoff; if x' = 0, B'B" is fair; 
and x' > 0, B'B" is losing, that is, has negative expected payoff. 

It is therefore possible to adjust the values of the aj and Pj in games B' and B" so 
that they are individually losing, but the combined game B'B" is now winning. This is 
the Parrondo effect. In the present example, the Parrondo effect occurs when 

(1 — a3)(l — 04) > aia2 (3.24) 

(l-/33)(l-/?4) > A/?2 (3.25) 

and 

il-h)il-U) <tit2. (3.26) 

The reader is referred to [15] for a detailed analysis of the values of the parameters 
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which lead to the Parrondo effect in such games. 

Restricting to the original work of Parrondo et al, a special case occurs when we 
consider one of the games in the randomized sequence to be of type A. That is, flipping 
a single biased coin which on the surface appears to have no history dependence. How- 
ever, note that such a game may be interpreted as a history dependent Parrondo game 
with a two stage history where the coin used in A is employed for every history. Call 
such a history dependent game A'. The transition matrix for A' takes the form 



A 



p 





p 





1 — p 





1 — p 








p 





p 





1 — p 





1 - 



(3.27) 



v 

Now, forming randomized sequences of games A' and B' is seen to agree with the 
forming of convex linear combinations mentioned above. In particular, as analyzed 
in [30] if games A' and B' are now sequenced randomly with equal probability, the 
Markov process for the randomized sequence is given with transition matrix containing 
the entries qj = ^(aj +p) and 1 — = ^[(1 — aj) + (l—p)] in the appropriate locations 
(recall that the probability of win for game A is p), and has stationary state 







P2 


1 




~ M 


P3 





( 



94(1 -gi) 
^4(1 -gi) 

\ (l-gi)(l-g2) I 



(3.28) 
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Denote this randomized sequence of games A' and B' by A'B'. The probability of gain 
in the game A'B' is 

M ^ ^^^^ (1 - gi) (2g4 + 1 - ^2) + g3g4 ^ 

As in the more general case of the game B'B", it is now possible to adjust the values 
of the parameters p and pj's in games A' and B' so that they are individually losing, but 
the combined game A'B' is now winning. This happens when 

l-p>p (3.30) 

(1 — a;3)(l — 0:4) > aia2 (3.31) 

and 

(l-g3)(i-g4) <gig2. (3.32) 

Parrondo et al show in [30] that when p = | — e, ai = ^ — e, 0:2 = 03 = | — e, 
a4 = — e, and e < j|g, the inequalities (3.30)-(3.32) are satisfied. This is Parrondo et 
al's original example of the Parrondo effect for history dependent Parrondo games. 

3.2 The FNA Quantization of Parrondo Games 

In [11], Flitney, Ng, and Abbott quantize the type A' Parrondo game by considering 
the action of an element of SU (2) on a qubit and interpret this as "flipping" a biased 
quantum coin. They consider history dependent games with (n — 1) stage histories, and 



45 



Chapter 3. Properly Quantizing History Dependent Parrondo Games 



in the language of the Bleiler formalism, quantize these games via a family of protocols. 
In every protocol, n qubits are required and the unitary operator representing the entire 
game is a 2" x 2" block diagonal matrix with the 2x2 blocks composed of arbitrary 
elements of SU (2). In the language of quantum logic circuits, this is a quantum multi- 
plexer [17]. The first (n — 1) qubits represent the history of the game via controls, as 
illustrated in Figure 3.1 for a two stage history dependent game similar to the game B' 
given in Table 3.2. Each protocol is defined as the action of the quantum multiplexer on 
the n qubits. 

The quantum multiplexer illustrated in Figure 3.1, where the elements Qi ... are 
elements of SU{2), operates as follows. When the basis of the state space (CP^)®'^ of 
three qubits is the computational basis 



the quantum multiplexer takes on the form of an 8 x 8 block diagonal matrix of the form 



B = {|000) , |001) , |010) , |011) , |100) , |101) , |110) , |111)} . 



Qi 



\ 



Q 







(3.33) 







\ J 
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03 



Figure 3.1: Part of the quantization protocol for the history dependent Parrondo game. The first 
two wires represent the history qubits. 

where each Qj E SU{2). That is 



aj -bj 



bj aj 



\ 
J 



(3.34) 



with aj, bj G C satisfying |aj|^ + \bjf = 1. 

For further description of the workings of the quantum multiplexer, the following 
convention, found in D. Meyer's original work [ ], will be used. Let a "win" or "gain" 
for a player be represented by the action "No Flip" which is the identity element of 
SU (2). For example, in Meyer's quantum penny flip game, the "quantum coin" is in the 
initial state of "Head" represented by |0) and a gain for the player using the quantum 
strategies occurs when the final orientation state of the coin is observed to be |0). This 
is contrast to the convention in FNA [. i] where |1) represents a gain. 

Now the first two qubits of an element of B represent a history of the classical game, 
with |0) representing gain (G) and the |1) representing loss (L). The blocks Qj act on 
the third qubit in the circuit under the control of the history represented by the binary 
configuration of the first two qubits. For example, if the first two qubits are in the joint 
state 1 00), the SU (2) action Qi is applied to the third qubit. Similarly, for the other three 
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basic initial joint states of the first two qubits. This models the historical dependence 
of the game by having the history (G, G) correspond to the initial joint state |00) of the 
first two qubits, the history (G, L) correspond to the initial joint state |01), the history 
{L, G) correspond to the initial joint state |10), and the history (L, L) correspond to the 
initial joint state |11). Thus, an appropriate action is taken for each history. 

Recall from section 3.1.2 that the evaluation of the behavior of the classical history 
dependent Parrondo game requires more than just the Markov process. The evaluation 
also requires the stable state and a payoff rule. Note that the results of applying the quan- 
tum multiplexer depends entirely on the initial state on which it acts. That is, different 
initial states result in differing final states. The payoff rule used by Abbott, Flitney, 
and Ng resembles that for the classical game in that the quantized versions are winning 
when the expectation greater than (gain capital), fair if the expectation is equal to 
(break even), and losing if the expectation is less than (lose capital). Further, as in the 
classical game this question is decided by examining the probability of gain versus the 
probability of loss. In particular, if the probability of gain is greater than ^, the quantum 
game is winning. 

3.2.1 Problems with the FNA Protocol 

The FNA quantization protocols for the history dependent game attempt to replace the 
classical biases of the coins in the game with arbitrary elements of SU{2) and the stable 
state of Markov process describing the dynamics of the game with certain initial states 
of the qubits on which a quantum multiplexer, composed of the arbitrary elements of 
SU (2), acts. The problems with the FNA quantization protocols are two-fold. First, the 
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attempted embedding of the classical history dependent game into the quantized game 
by replacing the biases of the classical coins with arbitrary SU{2) elements, turns out to 
be relational rather than functional. That is, Equations (3.33) and (3.34) together give 
a family of quantum multiplexers that the classical game maps into via the embedding. 
This relational mapping makes it impossible to recover the classical game by restricting 
the quantized game to the image of the embedding. Therefore, the FNA quantization of 
the history dependent Parrondo game is not proper. 

The second problem arises from the choice of initial state. No attempt is made 
to produce an analog of the stable state of a Markov process. Instead, the authors 
mention the obvious fact that different initial states will produce different results, and 
in particular consider two arbitrary initial states, one the maximally entangled state 

(|000) + |111)), the other the basic state |000). In the latter, the authors assert that 
the quantum game behaves like a classical game with fixed initial history (L, L), ac- 
cording to their convention in which |0) represents loss. Note that even if the this is not 
a proper quantization of any classical history dependent game as it fails to incorporate 
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|000) 



the other histories represented in the stable state. For 

(A 








and when acted upon by the quantum multiplexer in Equation (3.33) produces the output 

hi 






V / 

which makes the failure of the protocol to incorporate the other histories apparent. 

In the former, a similar situation occurs where only the histories |000) and |111) are 
incorporated. This protocol is also not proper as only the histories (L, L) and (G, G) are 
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non-trivially represented in the initial state. For 



V2 



(|000) + |111)) 



1 

71 



ll\ 








and when acted upon by the quantum multiplexer in Equation (3.33) produces the output 

( a. ^ 

hi 





\ ai J 



1 

71 



from which, again, the failure of the protocol to incorporate the other histories is appar- 
ent. 

Moreover, both quantization protocols fail to reproduce the Markovian dynamics 
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and the payoff function of the original game. 

Flitney et al also consider various "sequences" of the quantum games A' and B', 
where B' is played with three qubits and quantized using the maximally entangled initial 
state. These sequences are defined by compositions of the unitary operators defining the 
games. Indeed, these sequences now produce the results presented in [ 1 1]. These results 
are certainly novel and perhaps carry scientific significance; however, they fail to carry 
game-theoretic significance as, with respect to the classical Parrondo games, each arises 
from a quantization that is not proper. 

In light of the Bleiler formalism discussed in chapter 2, constructing proper quan- 
tizations of games is a fundamental problem for quantum theory of games. In the fol- 
lowing section, a proper quantization paradigm is developed for both history dependent 
Parrondo games and randomized sequences of such. 

3.3 Properly Quantizing History Dependent Parrondo Games 

Consider the history dependent game B' with only 2 histories. As in the FNA protocol, 
the quantization protocol for this game uses a three qubit quantum multiplexer with 
matrix representation 



Q = 



f 


Qi 











\ 







Q2 


















Qs 







\ 











Q4 


/ 



with each Qj G SU{2), together with an initial state. 

To reproduce the classical game, first embed the four classical coins that define the 
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game B' into blocks of the matrix Q corresponding to the appropriate history. The 
embedding is via superpositions of the embeddings of the classical actions of "No Flip" 
and "Flip" on the coins into SU (2) given either by 



(3.35) 



N = 




"1 






-n 




1° 








" J 



or by 



N* 



i 
i 



\ 
I 



(3.36) 



—IT] 

if] 

with r]^ = 1. Call the embeddings in equations (3.35) basic embeddings of type 1 
and the embedding in equations (3.36) basis embeddings of type 2. Choosing the basic 
embedding of type 1 embeds the j^^ coin into SU{2) as 



Qj = Vp]n+J{i-p,)f 



-y/i-Pjf] 

y^/i- PjV vpj 



(3.37) 



where pj is the probability of gain when the j^^ coin is played in the classical game 
B' given in Table 3.2. Note that the probabilities pj of gaining are associated with 
the classical action in line with Meyer's original convention from [ ] where |0) 
represents a gain. Hence, the elements of the subset 

W= (|000),|010),|100),|110)) 



of B all represent possible gaining outcomes in the game. The probabihty of gain in the 
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quantized game is therefore the sum of the coefficients of the elements of W that result 
from measurement. 

Next, set the initial state / equal to 



En 
.7 = 1 



7r„ 



V / 



(3.38) 



where the tTj are the probabilities with which the histories occur in the classical game, 
as computed from the stationary state of the Markovian process of section 3.1.2. The 
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quantum multiplexer Q acts on I to produce the final state 



S7 = l 



^a/(1 -P2)t^2 



\ 



(3.39) 



?7a/(1 -P3)7r3 

Measuring the state F/ in the observational basis and adding together the resulting coef- 
ficients of the elements of the set W' gives the probability of gain in the quantized game 
to be 

(3.40) 



QB' 



which is equal to the probability of gain in the classical game. 

This proper quantization paradigm is based on the philosophy first discussed at the 
end of chapter 2. That is, a proper quantization of a classical game Gs that depends on 
an initial state s requires that s be embedded into an initial state / on which the quantum 
multiplexer acts. Here, the initial state s = (tti, 7r2, tts, 774) G [0,1]^^ embeds as the 
initial state / G (CP^)®^ given in expression (3.38). The resulting game Gf' is the 
quantization of the classical game Gs by the protocol 6/ which maps the tuple {Q,I), 
withg = {Qi,Q2,Q3,Qa) e [SU {2)]''Uo Fj G (CPi)®^ gi^g^ in Equation (3.39). 
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Figure 3.2: Proper Quantization, using the embedding e, of the History Dependent Game via 
the quantization protocol Qj. 

Formally, 

ej:{Q,I)^Fi (3.42) 

By projecting on to the gaining basis W, one now gets a quantum superposition over the 
image ImG of the game G. Finally, quantum measurement produces ImC Call Proj 
the function that projects Fj on to W, and denote quantum measurement by Qmeas- Then 

Gf' = qmeas o Pvoj o 0, : (Q, /) ^ ImG (3.43) 

is a proper quantization of the payoff function of the normal form of classical history 
dependent game Gs given in Equations (3.15) and (3.16). Equation (3.43) can be ex- 
pressed by the commutative diagram of Figure 3.2, which the reader is urged to compare 
and contrast with Figure 2.4 in chapter 2. 

Note that by embedding s into /, the notion of randomization via probability distri- 
butions is generalized in the quantum game to the higher order notion of randomization 
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via quantum superpositions plus measurement. In particular, the probability distribution 
P = (Pi;P2,P3,P4) G [0, 1]^"^ that defines the Markov process associated with the game 
is replaced with the quantum multiplexer Q = iQi,Q2,Q3,Q'i) £ [S't/(2)]^^ associated 
with the quantized game, and the stable state s of the Markov process is replaced with 
an initial evaluative state / of the quantum multiplexer. 

3.4 Properly Quantizing Randomized Sequences of History Dependent Parrondo 



Recall from section 3.1.3 that randomized sequences of games B' and B" are analyzed 
via a Markov process with transition matrix equal to a real convex combination of the 
transition matrices of each game in which B' is played with probability r and B" with 
probability (1 — r). Moreover, such a sequence is considered to by an instance of a 
history dependent game denoted as B'B". 

Motivated by the discussion on proper quantization of the game Parrondo games 
B' and B" in section 3.3 above, let us now consider a higher order randomization in 
the form of a quantum superposition of the quantum multiplexers used in the proper 
quantization of the the games B' and B" with the goal of producing a proper quantization 
of the game 

As in section 3.3, associate the quantum multiplexer Q' = {Q'i,Q'2, Q's, Q'^) with the 
game B', where 



Games 




\ 



I 
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Next, associate the quantum multiplexer Q" — {Q'(, Q2, Q'^, Q'[) with the game B", 



where 



Q';^^jN* + ^il-Pj)F* = 

Now consider the quantum superposition 

E = j'Q' + i'Q" 
( 



iQ'i + i'Q'l 000 

i'Q2 + i"Q2 

iQs + YQs 

V yg^ + i'Q'i J 

of the quantum multiplexers Q' and Q" with 



(3.44) 



(3.45) 



(7')' + (7T = l, W\' = r, |7"|' = (l-r), 7 '7" - 7'V = (3.46) 



and 



iQ'+i'Q'' = 



/ 

(3.47) 
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Set the evaluative initial state in this case equal to 



r 



\ 







T. 



2 



1 







(3.48) 









V 







where the Tj are the probabilities that form the stationary state of the classical game 
B'B" given in Equation (3.19). The claim is that the quantum multiplexer S in Equa- 
tion (3.44) together with the evaluative initial state / in Equation (3.52) define a proper 
quantization of the classical game B'B" in which B' is played with probability r and 
and B" is played with probability (1 — r). 

To check the validity of this claim, compute the output of S for the evaluative initial 



59 



Chapter 3. Properly Quantizing History Dependent Parrondo Games 



state / in Equation (3.52): 



En 



v^3(7'v^3 + 7V^0 
The probabihty of gain produced upon measurement of this output is 



1 I 



(3.49) 



which simplifies to 



^ 4 

^ [l^'l^ + l7'f /^j + V^jPji (77" - 7^7') 



(3.50) 



Using the conditions set up in Equation (3.46), the previous expression further simplifies 
to give 

1 ^ 1 ^ 



R 



which is exactly that given in Equation (3.20) in section 3.4 for the classical game B'B" . 
Again, note that this proper quantization paradigm requires mapping of the initial 
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state of the classical game B'B", which is a probability distribution, into an initial state 
which the quantization protocol acts on, which is a higher order randomization in the 
form of a quantum superposition which measures appropriately with respect to the ob- 
servational basis. The image of the normal form of the quantum game in [0, 1] agrees 
precisely with p^^:^^ . Note that in this proper quantization of B B , not only is the 
initial state of the classical game replaced by a quantum superposition, but also a prob- 
abilistic combination of the transition matrices of the classical games is replaced with a 
quantum superposition of the quantum multiplexers associated with each classical game. 

3.4.1 A Special Case 

Recall from section 3.1.3 the classical analysis of the special case of the randomized se- 
quence of history dependent Parrondo games, with r = (1 — r) = ^, in which one of the 
games is A . The game A has the property that regardless of history, game A is always 
played. Such a sequence was considered to by an instance of a history dependent game 
denoted by A B' . In this section, a proper quantization of the randomized sequence is 
shown to follow as a special case of the proper quantization of the classical game B'B" 
developed in section 3.4 above. 

As before, associate the quantum multiplexer Q' = {Q[, Q'2, Q'^, Q'^), where 



with the game B'. Now, first embed the game A into SU(2) using basic embeddings of 
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type 2. That is, 



\ 



yfpi -y/l-p{ir]) 



\ 
I 



A 



The transition matrix for the game A' was given in Equation (3.27) and is reproduced 
here: 

^ p p 
l-p l-p 
p p 
y l-p l-p J 

The form of A suggests that the quantum multiplexer Q" = {A, A, A, A) should be 
associated with the game A'. Now let 7' = 7" = Equation (3.44) so that 



V2 



(A' + Q') 



1 

71 



[ A + Q\ ^ 

A + 

A + Q'3 

V A + Q'^ J 



(3.51) 
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with 



A + g' 



( 



\ 



(v^i ~Pj + VI - pi) V 



- iV^-piiv) + v^l - PjV) 

-- (v^i -Pj - VI - ^ 



With the evaluative initial state 



En 
,=1 Pi 



v 



/ 



/Pj- Vp^ 



J 





VP2 



VPs 


VP4 





(3.52) 



where the pj are the probabilities that form the stationary state of the classical game 
A B' given in Equation (3.28), the quantum multiplexer S in Equation (3.44) defines 
a proper quantization of the classical game AB' when both A and B' are played with 
equal probability. 

To see this, compute the output of S for the evaluative initial state / in Equation 
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(3.52): 



VPi (VI - Pi + VI -p-i) V 

Vp2(Vp^ + VP2) 
V^2 (VI -P2 + VI - pi) V 
2 E^i pj V^3 ( Vp^ + V^) 

VPs (VI ~P3 + VI -p^) V 

\ VP4 (VI -P4 + VI - p^) ^ y 

The probabihty of gain produced upon measurement is 



1 fP + Pj\ 1 A 



(3.53) 

which is exactly that given in Equation (3.29) in section 3.1.2 for the classical game 
A'B'. 



3.5 A Second Proper Quantization of the Randomized Sequence of History De- 
pendent Parrondo Games 



A second proper quantization of the sequence B'B" can be constructed in a manner 
similar to that used to construct the proper quantization for B' in section 3.3. Instead 
of forming a quantum superposition of the quantum multiplexers associated with each 
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game, first embed the classical coins used in the game B' B" into SU (2) as 

with 

=ra, + (l-r)/3, and 1 - = r(l - a,) + (1 - r)(l - 

and associate the quantum multiplexer Y = (^1,^2,^3,^) with the classical game 
B'B" . Set the initial state, as in section 3.4, equal to 
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where the r, are the probabilities that form the stationary state of the classical game 
B'B" given in Equation (3.19). The output state of this protocol is 



En 



A/r2(l -t2)ri 



\ V^4(l -U)r] J 



which, upon measurement produces the probability of gain 



(3.54) 



QB'B" _ 1 , 
Pgain ~ ■spn / , '^3 '^3 

^3=1 j=l 



which is exactly the probability of gain computed in Equation (3.29) of section 3. 1 .2 for 
the classical game AB'. 

Hence, there are two approaches, both motivated by different facets of the Bleiler 
formalism, used here to properly quantize random sequences of Parrondo games A and 
B' in which each game occurs with equal probability. One approach, discussed in sec- 
tion 3.3, generalizes the notion of randomization between the two games via probability 
distributions to randomization between games via quantum superpositions. The other 
approach, discussed above, embeds a probabilistic combination of the games into a 
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quantum multiplexer directly rather than via quantum superpositions of the protocols 
for each game. 

In the former approach, note that it was crucial that game A was embedded into 
SU (2) using basic embedding of type 2 as this allowed for the use of the broader arith- 
metical properties, namely factorization, of complex numbers to reproduce the classical 
result. In the latter on the other hand, basic embedding of type 1 sufficed. 

These two different approaches to quantizing history dependent Parrondo games 
raise interesting questions regarding the relationship between general quantum history 
dependent Parrondo games, which are quantum multiplexers with arbitrary SU (2) ele- 
ments forming the diagonal blocks, and the proper quantizations of the classical history 
dependent Parrondo games. For example, can a general quantum history dependent Par- 
rondo game always be factored into a sum of games which correspond to embedding of 
some classical history dependent Parrondo games? The reader is referred to the future 
directions section of chapter 6 where this subject is discussed in detail. 
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QUANTUM LOGIC SYNTHESIS BY DECOMPOSITION 



In Chapter 3, quantum multiplexers were used to properly quantize certain history de- 
pendent Parrondo games. In the following, quantum multiplexers will play a central role 
in synthesis of quantum logic circuits. 

Recent research in generalizing quantum computation from 2-valued qubits to d- 
valued qudits has shown practical advantages for scaling up a quantum computer. A 
further generalization leads to quantum computing with hybrid qudits where two or 
more qudits have different finite dimensions. Advantages of hybrid and d-valued gates 
(circuits) and their physical realizations have been studied in detail by Muthukrishnan 
and Stroud [ ], Daboul et. al [ ], and Bartlett et. al [2]. 

Recall from section 1.1 that the evolution of state space changes the state of the 
qudits under the action of a unitary matrix. Because evolution matrices are viewed as 
quantum logic gates in quantum computing, an essential idea from the theory of classical 
logic circuits carries over, namely, logic synthesis. One of the goals of logic synthesis 
is to express a given logic gate in terms of a universal set of quantum logic gates. Re- 
call from section 1.2.2 that sets of one and two qubit (even qudit) gates are universal. 
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Hence, the synthesis of a quantum logic gate requires that the corresponding matrix be 
decomposed to the level of unitary matrices acting on one or two qudits. Technological 
considerations for the implementation of one qudit gates might still require synthesis of 
these gates in terms of simpler one qudit rotation gates and two qudit controlled rotation 
gates. For 2-valued quantum computing, this is easily accomplished by the well known 
Euler angle parameterization of a 2 x 2 special unitary matrix (since a unitary matrix 
is equivalent to a special unitary matrix up to a complex multiple). For higher valued 
quantum computing, Tilma et al's work in [35] shows that a one qudit gate can be syn- 
thesized in terms of an Euler angle parametrization similar to the one available for 2 x 2 
special unitary matrices. 

If the quantum system consists of multiple qudits, then a gate may be synthesized by 
matrix decomposition techniques such as QR factorization and the cosine- sine Decom- 
position (CSD). Both the acronym CSD and the term CS decomposition will be used to 
refer to the cosine-sine decomposition from now on. The CSD is used by Mottonen et. al 
[22] and Shende et. al [31] to iteratively synthesize multi-qubit quantum circuits. Khan 
and Perkowski [16] use the CSD to develop an iterative synthesis method for 3- valued 
quantum logic circuits acting on n qudits. Bullock et. al present a synthesis method 
for n qudit quantum logic gates using a variation of the QR matrix factorization in [7] 
In [17], Khan and Perkowski give a CSD based method for synthesis of n qudit hybrid 
and (i-valued quantum logic gates. This chapter reviews the work of these authors on 
quantum logic synthesis techniques based on the CS decomposition. 



69 



Chapter 4. Quantum Logic Synthesis by Decomposition 



4.1 The Cosine-Sine Decomposition (CSD) 

Let the unitary matrix W G C'"^™ be partitioned in 2 x 2 block form as 



W 



m — r 



r 


m — r 






W21 


W22 



(4.1) 



with 2r < m. Then there exist r x r unitary matrices U and X, r x r real diagonal 
matrices C and S, and (m — r) x (m — r) unitary matrices and Y such that 



u 
V 



C -S 
S C 
J, 



m— 2r 



V 



Y 



(4.2) 



The matrices C and 5 are the so-called cosine-sine matrices and are of the form C 
= diag(cos 61, cos 62, ... , cos 6r), S = diag(sin 61, sin 62, ... , sin 6r) such that sin^ 6i + 
cos^ = 1 for some 9i, I < i < r [ ]. Algorithms for computing the CSD and 
the angles 9i are given in [4, 33]. The CSD is essentially the well known singular value 
decomposition of a unitary matrix implemented at the block matrix level [28]. Appendix 
B gives a review of the CS decomposition. 

The reader is advised that in the narrative that follows quantum logic gates, circuits 
and the corresponding unitary matrices will not be distinguished. 
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4.2 Synthesis of 2- valued (binary) Quantum Logic Circuits 



As the authors of [16, 22, 31, 36] show, CSD gives a recursive method for synthesizing 
2-valued and 3-valued n qudit quantum logic gates. In the 2-valued case the CSD of a 
2" X 2" unitary matrix W reduces to the form 



W 



U 
V 



\ 



J 



( 



\ 



c -s 
s c 



v 



Y 



(4.3) 



with each block matrix in the decomposition of size 2"~^ x 2"~^ 

A quantum multiplexer is a quantum logic gate acting on n qubits of which one is 
designated as the control qubit. If the control qubit of a quantum multiplexer is the 
lowest order qubit, that is, the first qubit in the joint state of n qubits, the multiplexer 
matrix is block diagonal. Note that the lowest order qubit is represented as the top most 
qubit in circuit diagrams. Thus, in terms of synthesis, the block diagonal matrices in 
Equation (4.3) are quantum multiplexers [31]. Now, depending on whether the control 
qubit carries 1 0) or 1 1 ) , the gate then performs either the top left block or the bottom right 
block of the n x n block diagonal matrix on the remaining (n — 1) qubits, respectively. 
A circuit diagram for a n qubit quantum multiplexer with the lowest order control qubit 
is given in Figure 4.1 where the black circle represents control via the basis state |1). 

Such a quantum multiplexer is expressed as 



|ai) (g) 



Uo 
Ui 



\ 



(|a2) ® ■ ■ ■ ® |a„)) 



(4.4) 
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+1 



+1 



M 



Figure 4.1: 2-valued Quantum Multiplexer M controlling the lower {n — 1) qubits by the top 
qubit. The slash symbol (/) represents (n — 1) qubits on the second wire. The gates labeled +1 
are shifters (inverters in 2-valued logic), increasing the value of the qubit by 1 mod 2 thereby 
allowing for control by the highest qubit value. Depending on the value of the top qubit, one of 
Ut is applied to the lower qubits for t G {0, 1}. 

where |aj) is the i-th qubit in the circuit, and both block matrices Uq and Ui are of size 
2n-i ^ 2""^. Depending on whether |ai) = |0) or |ai) = the expression (4.4) 
reduces to 

|0)®?7o(|a2)® ks)®---® kn)) (4.5) 



or 



|l)®f/i(|a2)® las)®---® |a„)) 



(4.6) 



respectively. 

A uniformly {n — 1) -controlled Ry rotation gate Ry is composed of a sequence of 
(n — l)-fold controlled gates Ry\ all acting on the lowest order qubit, where 



R^' 



cos y,- — sm 6',- 



sin Qi cos Qi 



(4.7) 



The cosine-sine matrix in Equation (4.3) is realized as a uniformly (n — l)-controlled 
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R'; 

























o 



o- 



-o- 



-o- 



Figure 4.2: A uniformly (n — 1) -controlled i?^^ rotation for 2-valued quantum logic. The o 
control turns on for control value |0) and the • control turns on for control value It requires 
2"~^ one qubit controlled gates Ry' to implement a uniformly (n — l)-controlled Ry rotation. 

Ry rotation gate, a variation of a quantum multiplexer, as shown in Figure 4.2. The 
control selecting the angle 9i in the gate Ry^ depends on which of the (n — 1) basis state 
configurations the control qubits are in at that particular stage in the circuit. In Figure 
4.2, the white circle represents control via the basis state |0). The i-th [n — l)-controlled 
gate Ry^ may be expressed as 



COS 6i — sm Uj 



sin di cos Oj 



ki) ® (|ai) ® ■■■ ® |a„)) 



(4.8) 



with 6i taking on values from the set {^^o, ^i, • • • , ^^2"-i-i} depending on the specific 
configuration of (|a2) ® ■ ■ ■ ® |a„)), resulting in a specific Ry^ for each i. 

As an example, consider the 3 qubit uniformly 2-controlled Ry gate controlhng the 
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-O 



+1 



Figure 4.3: A control by input value (mod 2) realized in terms of control by the highest value 
1 (mod 2). 













R"; 




























1 



-o 



o — 

Figure 4.4: A uniformly 2-controlled Ry rotation in 2-valued logic: the lower two qubits are 
the control qubits and the top bit is the target bit. 

top qubit from Figure 4.4. Then the action of R^^ on the circuit is 



COS Ui — sm 6',- 



v 



sin 9j cos 9i 



kl) ® (1^2) ® Wi)) 



(4.9) 



with 9i E {9o, 9i, 92, 9^}. As \a2) ® jaa) takes on the values from the set 

{|0) (g) |0) , |0) (g) |1) , |1) (g) |0) , |1) (g |1)} in order, the expression in (4.9) reduces to 

the following 4 expressions respectively. 



cos c^o — sm 9o 
^ sm9o cos 6*0 y 



|ai)®(|0)® |0)) 



(4.10) 



cos 9i — sm 9i 
sin 6*1 cos 6*1 



|ai)®(|0)® |1)) 



(4.11) 



74 



Chapter 4. Quantum Logic Synthesis by Decomposition 



cos ^2 ~ sm 



v 



sin Qi cos Qi 



(4.12) 



V 



|ai)®(|l)® |1)) 



(4.13) 



cos 6*3 — sin Q-i 
sin ^3 cos 6'3 

Observe that by iterating the CSD and factoring the result each time results in a 
quantum circuit consisting of variations of the quantum multiplexer. 



4.3 CSD Synthesis of 3-valued (ternary) Quantum Logic Circuits 



In the 3-valued case, two applications of the CSD are needed to decompose a 3"^ x 3" 



unitary matrix W to the point where every block in the decomposition has size 3 



n-l 



jn— 1 



[16]. Choose the parameters m and r given in Equation (4.1) as m = 3" and 



r = 3"-\ so that m - r = 3" - 3""^ = 3"-^(3 - 1) = 3""^ ■ 2. The CS decomposition 
of W will now take the form in Equation (4.2), with the matrix blocks U and X of size 
3"^^ X 3"^^ and blocks V and Y of size 3""^ ■ 2 x 3"^^ • 2. Repeating the partitioning 
process for the blocks V and Y with m = 3"^^ • 2 and r = 3"^^, and decomposing 
them with CSD followed by some matrix factoring will give rise to a decomposition of 
W involving unitary blocks each of size 3""^ as follows. 



W = ABC 



( C -5 ^ 



s c ^ 

0/ 



DBF 



(4.14) 
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+2 



M 



+2 



+2 



u. 



Figure 4.5: 3-valued Quantum Multiplexer M controlling the lower (n — 1) qutrits via the top 
qutrit. The slash symbol (/) represents (n — 1) qutrits on the second wire. The gates labeled 
+2 are shift gates, increasing the value of the qutrit by 2 mod 3, and the control turns on for 
input 1 2). Depending on the value of the top qutrit, one of Ut is applied to the lower qutrits for 
tG {0,1,2}. 



-O 



+2 ^ -2 



Figure 4.6: A control by the value (mod 3) realized in terms of control by the highest value 2 
(mod 3). 



with 



A 



D 



( 



\ 



























Yi 

























B 



E 



/ 

Ci -Si 

V 5i Ci y 

/ ^ 

C2 -S2 
S2 C2 



\ 



J 



c 



^10 ^ 

Zi 

\0 Z2 J 

(4.15) 

^10 ^ 



v 



Wi 
W2 



(4.16) 

We realize the block diagonal matrices A,C,D and F in (4.15) and (4.16) as 3- 
valued quantum multiplexers acting on n qutrits of which the lowest order qutrit (top 
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Figure 4.7: A control by the value 1 (mod 3) realized in terms of control by the highest value 2 
(mod 3). 



-o- 



R 



Figure4.8: Auniformly (n — 1) -controlled rotation. The lower (n—1) qutrits are the control 
qutrits. The controls o, •, and o turn on for inputs |0), |1), and |2) respectively. It requires 3"^^ 
one qutrit controlled gates to implement a uniformly (n — 1) -controlled Rx or Rz rotation. 

most in a circuit diagram) is designated as the control qutrit. Depending on which of the 
values |0), |1), or |2) the control qutrit carries, the gate then performs either the top left 
block, the middle block, or the bottom right block respectively on the remaining n — 1 
qutrits. Figure 4.5 gives the layout for a n qutrit quantum multiplexer realized in terms 
of Muthukrishnan-Stroud (MS) gates. The MS gate is a rf-valued generalization of a 
controlled gate from 2-valued quantum logic, and allows for control of one qudit by the 
other via the highest value of a (i-valued quantum system, which in the 3-valued case is 
2 [23]. 

The cosine-sine matrices are realized as the uniformly [n — 1) -controlled and Rz 
rotations in M?. Similar to the 2-valued case, each R^ and R^ rotation is composed of a 
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sequence of (n — l)-fold controlled gates i?^.' or Rf \ where 



1 

COS 9i — sin 9i 
sin 6i cos 6i 



J 



\ 



cos (pi 


— sm 


sin (pi 


cos q 









(4.17) 



Each i?^' or i?^' operator is applied to the top most qutrit, with the value of the angles 
6i and 0i determined by the (n — 1) basis state configurations of the control qutrits. A 
uniformly controlled Rx gate is shown in Figure 4.8. Figures 4.6 and 4.7 explain the 
method to create controls of maximum value. Notet that the value of the control qubit is 
always restored in Figures 4.6 and 4.7. 



4.4 Synthesis of Hybrid and rf- valued Quantum Logic Circuits 

It is evident from the 2 and 3-valued cases above that the CSD method of synthesis is of 
a general nature and can be extended to synthesis of (i-valued gates acting on n qudits. 
In fact, it can be generalized for synthesis of hybrid n qudit gates. We propose that a 
{did2 ■ ■ ■ dn) X {did2 ■ ■ ■ dn) block diagonal unitary matrix be regarded as a quantum 
multiplexer for an n qudit hybrid quantum state space 7i = Hd^ ® 'Hd2 ® ■ ■ ■ ® 'Hd^, 
where Hd, is the state space of the i qudit. 

Moreover, consider a cosine-sine matrix of size {did2 ■ ■ ■ dn) x {did2 . . . dn) of the 
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form 



Jp ^ 



(4.18) 



C -S 
S C 
V /, / 

with Ip and Jg both some appropriate sized identity matrices, C = diag(cos 9i, 
cos 02, ... , COS 9t) and S = diag(sin ^^i, sin 6*2, ... , sin 9t) such that sin^ 9i + cos^ 9i = 1 
for some 9i with 1 < i < t, and p + q + 2t = {did2 . . . dn). We regard this matrix as 
a uniformly controlled Givens rotation matrix, a generalization of the Ry, R^, and Rz 
rotations of the 2 and 3-valued cases. A Givens rotation matrix has the general form 



... cos 9 ... — sin 9 ... 







sin^ 







cos 6* 







(4.19) 



where the cosine and sine values reside in the intersection of the i-th and j-th rows 
and columns, and all other diagonal entries are 1 [ ]. Hence, a Givens rotation matrix 
corresponds to a rotation by some angle 9 in the ij-th hyperplane. 

Based on the preceding discussion, we give in Theorem 4.4.1 below an iterative 
CSD method for synthesizing a n qudit hybrid quantum circuit by decomposing the 
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corresponding unitary matrix of size {did2 . . . dn) x {did2 ■ ■ ■ dn) in terms of quantum 
multiplexers and uniformly controlled Givens rotations. As a consequence of Theorem 
4.4.1, we give in corollary 4.4.1 a CSD synthesis of a quantum quantum logic circuit 
with corresponding unitary matrix of size rf" x c/". The synthesis methods given above 
for 2-valued and 3-valued circuits may then be treated as special cases of the former. 

4.4.1 Hybrid Quantum Logic Circuits 

Consider a hybrid quantum state space of a n qudits, Ti = Hdj ® Hd^ ® • • ■ ® 'Hd„, 
where each qudit may be of distinct d- valued dimension di, 1 < i < n. Since a qudit in 
His a column vector of length did2 . . . dn, a quantum logic gate acting on such a vector 
is a {did2 ■ ■ ■ d]\f) x {did2 . . . dn) unitary matrix W. We will decompose W, using CSD 
iteratively, from the level of n qudits to (n — 1) qudits in terms of quantum multiplexers 
and uniformly controlled Givens rotations. However, since the rf-valued dimension may 
be different for each qudit, the block matrices resulting from the CS decomposition may 
not be of the form d^~^ x d'^~^ for some d. Therefore, we proceed by choosing one of 
the qudits, of dimension rfj, to be the control qudit and order of the basis of H in 
such a way that is the highest order qudit. We will decompose W with respect to 
Q. so that the resulting quantum multiplexers are controlled by c^- and the uniformly 
controlled Givens rotations control c^^ via the remaining (n — 1) qudits. We give the 
synthesis method in Theorem 4.4.1 below. 

Theorem 4.4.1. Let be an M x M unitary matrix, with M = did2 ■ ■ - dn, acting as 
a quantum logic gate on a quantum hybrid state space H = Hd^ ® Hd^ ® ■ ■ ■ ® Hd„ of 
n qudits. Then W can be synthesized with respect to a control qudit Cd^ of dimension 
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di, having the highest order in H, iteratively from level n to level — 1) in terms of 
quantum multiplexers and uniformly controlled Givens rotations. 

Proof. Step 1. At level n, identify a control qudit of dimension d^. Reorder the basis 
of H so that Crf. is the highest order qudit and the new state space isomorphic to H is 

If we choose values for the CSD parameters m and r as m = {did2 ■ ■ ■ dn) and 
r = ((^1^2 • • • di_idi+i . . . dn), thenm — r = di. . . di_idi^i . . . dn{di — l). Decomposing 
W by CSD, we get the form in (4.2) with the matrix blocks U and X of size r x r and 
blocks V and Y of size (m — r) x (m — r). Should m — r not have the factor (di — 1), 
we would achieve the desired decomposition of W from level of n qudits to the level of 
(n — 1) qudits in terms of block matrices of size r x r. The task therefore is to divide 
out the factor (di — l) from m — r by an iterative lateral decomposition described below, 
that uses the CSD to cancel [d^ — 1) from m — r at each iteration level leaving only 
blocks of size r x r. 

For step 2 of the proof below, we will say that a matrix with k rows and k columns 
has size k instead oik x k. 

Step 2. Iterative Lateral Decomposition: For the unitary matrix W of size M, we 
define the j-th lateral decomposition of W as the CS decomposition of all block matrices 
of size other than r that result from the (j — l)-st lateral decomposition of W: 

For 0<j<{d^- 2), set 
mo = {did2 ...dn) 
fo = {did2 . . . di^idij^i . . . dn) 
If] = 
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Apply CSD to W 
Else set 

rrij = mo- j - To 
Tj = ro 

"b' ^ = "^0 - (i + l)'^o 

= {dxd2 . . . di-idi+i ...dn) [di - {j + 1)] 

ruj - 2rj = mo - {j + 2)ro 

= {did2 . . . di_idi+i ...dn) [di - {j + 2)] 

Apply CSD to matrix blocks of size other than rofrom step j — 1 
End If 
End For. 

When j — 0, we call the resulting 0-th lateral decomposition the global decomposi- 
tion. Note that if di — 2, then the algorithm for the lateral decomposition stops after the 
global decomposition. This suggests that whenever feasible, the control system in the 
quantum circuit should be 2- valued so as to reduce the number of iterations . Below we 
give a matrix description of the algorithm. 

For j = 0, the 0-th lateral decomposition of W will just be the CS decomposition of 

W. 

W = A^Bj^^D^^^ (4.20) 
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where 



A 



(0) 



(0) 



a 



(0) 



(0) 



(0) 



Vr 



(0) 



s: 



(0) 



(0) 






^mo-2ro y 



(0) 



X 



(0) 



Yr 



(0) 



with ?7q°\ Xq°\ C'^\ and 5*0° ■* all of the desired size ro, while Vq^^ and Y'q'^^'' are of size 
mo — To- The superscripts label the iteration step, in this case j = 0. The subscript is 
used to distinguish between the various matrix blocks U, V, X, Y, C, S, that occur at the 
various levels of iteration. The 0-th lateral decomposition in the form from Equation 
(4.20) is called the global decomposition of W. 

For j = 1, we perform lateral decomposition on the blocks V^'^'^ and Y^'^^ of the 
block matrices A^^^ and D^^^ respectively, the only blocks of size other than tq resulting 
from the 0-th lateral decomposition given in (4.20). In both cases, set mi = mo — ro 
and ri = ro so that mi — ri = mo — 2ro. For V^'-^-' this gives the decomposition 



A 



(0) 



\ 



v 



(1) 



V, 



(1) 







(1) 



v 







(1) 










^mo-3ro J 



X, 



(1) 



(1) 



(4.21) 



with f/(S^\ Xi^\ Cj^) and Sj^^ all of size ro, and vj^^ and Y^^^^ of size mi - ri. All three 
matrices residing in the lower block diagonal of the matrix (4.21) are the same size. 
Therefore, by introducing identity matrices of size ro and factoring out at the matrix 
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block level, A^^^ will be updated to 



/1q — /1q i3Q 



(4.22) 



where 



A, 



(1) 














^0 











y{l) 
^0 



^(1) 











































I mo- 



3ro / 



(1) 



xi^) 












(1) 



/ 

A similar lateral decomposition of the block Y^^^ will update D^^^ in (4.20) to 



(0) 



(4.23) 



where 



(1) 




























( 



(1) 







c« 
cj^) 

\ 
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-3ro / 
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(1) 



xf) 










For iteration j ^ 0, perform lateral decomposition on the total 2^ blocks V^'' ^\ Y^f ^\ 



where {] < k < 2^^ 



-1) 



1, that occur in the global decomposition at the end of iteration 



(j — 1). For each V^f ^\ ^\ set Vj = tq, rrij = rrij^i — rj^i = ttlq — jvQ. For each 
V^^~^\ the lateral decomposition at level j will give the following 



lO-i) 



/ jjU) 



( c^o) 



\ 







(i) 



qU) 






Imo~(j+2)ro ) 



\ 



Yt, 



U) 



(4.24) 



where the A'^-'^^^ is the block diagonal matrix of size of j ■ tq arising from the lateral 
decomposition in the previous j steps. The blocks ujj\ xl^,\ c\^) and S\^) are all of 
size ro, < A;' < 2-' — 1. The blocks VjiP and yI:P are of size rrij — r^. The three 
matrices residing in the lower block diagonal of the matrix (4.24) are all of same size. 
Therefore, by introducing identity matrices of size j ■ ro and factoring out at the block 
level, A^i^^^ will be updated to 



/lO-i) _ aU) tj(J) n(j') 

/ii. — /ii,/ -D,, ly,, 
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+(4-1) Upy^ +(4-1) ypy- -A +(4-1) 



u. 



M 



Figure 4.9: An n qudit hybrid quantum multiplexer, here realized in terms of Muthukrishnan- 
Stroud (d-valued controlled) gates. The top qudit has dimension dj and controls the remaining 
(n — 1) qudits of possibly distinct dimensions which are represented here by the symbol (/). The 
control turns on for input value \di — 1) mod di of the controlling signal coming from the top 
qudit.The gates +{di — 1) shift the values of control qudit by (cZj — 1) mod di. 



where 



A 



0-1) 













KV 



(j) 
















~'^k' 

























Imo-{j+2)ro ) 



^k' 
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^k' 



\ 



J 



For the next iteration, set k — k' and iterate. Upon completion of the lateral decomposi- 
tion, repeat steps 1 and 2 for the synthesis of the circuit for the remaining (n — 1) qudits, 
with the restriction that each gate in the remaining circuit be decomposed with respect 
to the same control qudit identified in step 1. □ 
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Figure 4.10: A hybrid uniformly (n — 1) -controlled Givens rotation. The lower (n — 1) qudits 
of dimensions d2, d^, ■ ■ ■ , di-i, dj+i, . . . , d„, respectively, are the control qudits, and the top is 
the target qudit of dimension di. The control gate d] turns on whenever the control qudit of 
dimension di takes on the value k (mod) di. 



Since the basis for Ti was reordered in the beginning so that the control qudit was of 
the highest order, the block diagonal matrices with all blocks of size tq x tq are inter- 
preted as quantum multiplexers and the cosine-sine matrices are interpreted as uniformly 
controlled Givens rotations. In Figures 4.9 and 4.10, we present the circuit diagrams of 
a hybrid quantum multiplexer and a uniformly controlled Givens rotation, respectively. 
A uniformly controlled Givens rotation matrix on n qudits can be realized as the com- 
position of various {n — l)-fold controlled Givens rotation matrices, G^(*ijy acting on 
the top most qudit of the circuit with the angle of rotation depending on the basis state 
configuration, in their respective dimensions, of the lower [n — 1) qudits. 
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4.4.2 (i- valued Quantum Logic Circuits 

Given the hybrid n qudit synthesis, the case of d- valued synthesis becomes a special case 
of the former since by setting all di = d, the state space H = Hd^ Cg> 'Hd2 ® ■ ■ ■ ® Ti-d^ 
reduces to the state space Hf"'. Unitary operators acting on the states in iff" are unitary 
matrices of size c?" x rf". We give the following result for rf-valued synthesis. 

Corollary 4.4.1. A d- valued n qudit quantum logic gate can be synthesized in terms of 
quantum multiplexers and uniformly controlled Givens rotations. 

Proof: Since all the qudits are of the same dimension, there is no need to choose a con- 
trol qudit. In the proof of Theorem 4.4.1, set di = d for all i. Then M = did2 . . . dn = 
d". For iteration j = of the lateral decomposition, set mo = rf", tq = so 
that mo — ro = d"'~^{d — 1). For < j < (d — 2), set rj = tq = d^~^, and 
rrij = rrij^i — rj_i = d^^^^d — (j + 1)). 

For the rf-valued case, we note that there are a total of d"'~^{2'^^^ — 1) one qudit 
Givens rotations in the circuit at the {n — 1) level, each arising from the ^j-lg^'' 2* = 
2^^^ — 1 uniformly controlled Givens rotations in the CS decomposition of an n qudit 
gate. Moreover, in each uniformly controlled Givens rotation, there are (n — 
control symbols of which {n — correspond to control by the highest value of 

d — I. The latter controls do not require shift gates around them to increase the value 
of the signal qudit to d — 1. Hence, there are (n — l)d"-^^ — [n — = [n — 

— control symbols that correspond to control by values other than d — 1 
and therefore need two shift gates (fig. 11) around them. This gives the total number 
of one qudit shift gates in each uniformly controlled rotation to be 2{n — — 
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whereby the total number of one qudit shifts and Givens rotations in the circuit 
at the {n - 1) level is 2{n - - d"-2)(2<^-i - 1) + ^"-^(2'^-^ - 1) = (S'^"^ - 

1) [2(n - - + 

There are 2*^"^ quantum multiplexers in the decomposition, each consisting of a total 
of 2d shift and controlled gates. Hence, there are a total of d • 2*^ one qudit and controlled 
gates in the (n — 1) level circuit. This gives a total, worst case, one qudit and controlled 
gate count in the circuit at level (n-1) to be (2'^-^-l) [2{n - - d"-^) + 

d-2'^. 
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A QUATERNIONIC CO-ORDINATIZATION OF BINARY QUANTUM 

COMPUTATION 



A quaternionic coordinatization of the players' quantum strategies in certain quantized 
games by Landsburg in [ ] gives him a computational framework for classifying poten- 
tial Nash equilibria in these games. This idea led Ahmed, Bleiler and Khan [ ] to con- 
struct a parallel coordinatization using octonions for another class of quantized games, 
giving the authors a computational framework for classifying potential Nash equilibria 
in these games. Motivated by these result, this chapter proposes a quaternionic coor- 
dinatization of binary quantum computation by putting quaternionic coordinates on the 
Lie group SU (2) of quantum logic gates acting on one qubit and on the projective com- 
plex state space CP^ of one qubit, with the eventual goal of providing an enhanced 
computational capability for circuit analysis. 

In general, one qubit quantum logic gates are unitary matrices with determinant 1 
or —1. However, a 2 x 2 unitary matrix is equivalent to a special unitary matrix up to a 
factor of i. That is, if 



U 



a 
\^ c d J 
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is a unitary matrix with determinant ad — cb — —1, then 



ia ib 



\ 



= I 



\ 



ic id 



J 



where t/' has determinant {ia){id) — {ic){ib) = —ad+cb = —l{ad — cb) = (— 1) = 
1 and is therefore special unitary. The factor i is regarded a unitary phase in any resulting 
calculations. For the remainder of this chapter, all instances of a unitary matrix U with 
determinant —1 will be replaced with its equivalent special unitary matrix U' e SU (2). 

Now by identifying both SU (2), the set of one qubit quantum logic gates, and CP^, 
the state space of a qubit, with unit quaternions Sp{l), we develop a quaternionic co- 
ordinatization of binary quantum computation. In this chapter, we will use the notation 
1 and J for the unit quaternions 1 and j respectively to emphasize their roles as control 
signals in the context of quantum computing. 

5.1 Identifying SU(2) with Sp(l) 

The Lie group 5'p(l) of unit quaternions can be considered as 

Sp{l) = {u = uol + MiJ : \uf = \uof + \u^f = + + (u'^^ + = l} • 
The Lie group SU (2) of 2 x 2 special unitary matrices is 



SU{2) 




\ 



J 



: a,/? e C and |q;| V |/?|^ = 1 
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The special unitary requirement suggests a strong connection between SU (2) and Sp{l). 
Indeed, we can set up a one to one correspondence between SU (2) and Sp{l) as follows. 
Consider EI as under the identification 



at + 131 



( \ 

a 



and let yil + y2^ eM and zil + e Sp{l). Recall that al = la for all a e C and 
form the product 

{yil + y2^){zit + 2^2 J) = yizil + y2Jzi + yiZ2^ + y2^Z2l 

= {yizi - 1/2^2)1 + {y2Zi + 1/1^2) J- 



Write this result as an element of via the identification as 



But 



V 



V 



yizi - y2Z2 
y2Zi + yiZ2 



yizi - y2Z2 
7/2^1 + yiZ2 



( 



Z\ —Z2 

Z2 Zl 



\ ( \ 

yi 



J 



yy2j 



so the result of the quatemionic product, as an element of C^, is in the image of the 
special unitary transformation 

/ - \ 

Zl —Z2 

Z2 Zl 
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acting on C^, establishing the following identification of SU (2) and Sp{l) 



( _ \ 



\ 



Z2 Zl 



Zi\l + Z2^. 



J 



In other words, right multiplication by a unit quaternion in EI corresponds to the action 
on the left of the corresponding special unitary matrix on C^. 

In fact, it is possible to make the quatemionic product compatible with the left action 
of a SU (2) element on C^. That is, the left action of a linear transformation on can 
be made to correspond to multiplication on the left by a unit quaternion in H by writing 
quaternions with scalars on the right. For then, we get 

(Izi + Iz2)(tyi + 1^2) = l-^iyi + ^Z2yi + zily2 + 1^2 Jy2 

= l{ziyi - ^21/2) + ^{z2yi + ziy^) 

which corresponds to 



V 



Ziyi - Z2y2 
Z2yi + ziy^ 



( 



Zl -Z2 

Z2 Zl 



\ ( \ 

yi 



J 



with the identification of SU (2) and Sp{l) given by 



/ _ \ 

Zl —Z2 

Z2 Zl 



IZi + IZ2. 
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Either of these two identifications of SU (2) with 5*^(1) introduces quaternionic co- 
ordinates on SU (2). We choose the latter due its salient property of keeping the quater- 
nionic product compatible with the left action of SU{2) on C^. In other words, we 
consider the quaternions as a right complex vector space. 

It is an easy check that this identification preserves multiplication in SU (2). If 



v 



a —(3 
(3 a 



\ 



( 



5 —7 
7 (5 



G SU{2), 



then their product in SU{2) results in 



a —[3 
P a 



—7 



7 6 



I 



^ a5 — iS'y —0:7 — /35^ 



\ 



(35 + 07 — /57 + a5 



which is identified with the unit quaternion 



l(a(5-/37) + J(/35 + a7), 



(5.1) 



while identifying the SU{2) elements with unit quaternions ^zr^f results in the quater- 
nionic product 

(la + J/5) (1(5 + J7) = l{a5- '^^) + J{p5 + 07) 

the result of which is exactly the quaternion in (5.1). In fact, this identification sets up a 
Lie group isomorphism between 5^(1) and SU (2). 
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5.2 Identifying Sp{l) with CP^ 

Observe that the Bloch sphere 

cpi = (c^ - {0}) /c* ^ syu{i) 

where C* = M+ x and 



X 


■1 




-\ 


( \ 

X 


[y j 











for x, y e C, both not equal to zero, and the scalar A e and is called phase. Note 
that we scalar multiply elements of CP^ on the right rather than the left, a conven- 
tion that is necessary for differentiating between scalar multiplication and the action of 
SU (2) on CP^ under the identifications. 

The Hopf map H : ^ CP^ is defined here as 



H : 



X 



\y ) 



yx 



-1 



with ^ considered to be the number |. On the Bloch sphere, the pure states are repre- 
sented by J and ^ corresponding to the vectors 



(A 



Vv 
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respectively. In general, ^ = jSa ^ corresponds to the vector 



a 



and up to unitary phase 



1 ^ 


\ 




\ 

a 


^ f3X 


J 


v 





We identify this element of the Bloch sphere with a unit quaternion representing its orbit 
in S'\ That is. 



X 



\y J 



where 



1, 



/ ■ \ 



I, 



Ix + ,Jy 



(5.2) 



J, 



^ ^ 



K 



(5.3) 



is the identification of the basis of (hence CP^) with the basis of EI as complex vector 
spaces. The identifications in equations (5.2) and (5.3) induce a product between ele- 
ments of SU (2) and elements of CP^ via quatemionic multiplication that is consistent 
with the left action of an appropriate SU (2) element on the elements of CP^. That is, 
for 

a a —p 

A= ,A= gCP\ eSU{2) 



P a 



J 
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quatemionic multiplication gives the product ★ between A and A as follows. 



A 



/ \ / A^ 

a 



+ J/3)(15 + J7) 



V 



aS — (3 J 
(35 + a7 



v 



a —(5 
(5 a 



\1 J 



5.2.1 Action of U{1) on CP^ 



Note that the unit complex numbers U (1) can be embedded into SU (2) via 



a 



a 



a 



and in this form act on CP^ as linear transformation instead of scalar multiplication. 
Our identifications respect this fact, as the following example shows. 



v 



Q! \ 5 

^ ^ ) \ l ) J 



t{a5) + J(a7) = at5 + q;J7 
= + J7). 



Note that even though in the expression a{t5-\- ck J7) the complex number a appears on 
the left, it does not represent scalar multiplication because of our convention that scalars 
multiply on the right. In fact, it's occurrence on the left of the quaternion 15 + aJ^ tells 
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us that it represents the action of as a linear transformation under the embedding 
in SU (2). 
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FUTURE DIRECTIONS 



The proper quantization protocols developed in chapter 3 for history dependent Par- 
rondo games using certain quantum multiplexers lend a game theoretic perspective to 
the study of quantum logic circuits via quantum multiplexers. Indeed, the notion of the 
Parrondo effect is now attached to quantum circuits, and it is natural to raise the follow- 
ing question: 1) can a genuine "quantum Parrondo effect" be characterized in quantum 
circuits through this game theoretic perspective? 

Moreover, to date there is no agreement in the literature on exactly what a quantum 
Markov process is. One difficulty lies in coming up with an appropriate definition of 
the "quantum" stable state. Our quantizations of history dependent Parrondo games are 
essentially specific quantized Markov processes involving specific elements of the Lie 
group SU (2) and with stable states chosen game-theoretically. A more general set up 
is possible in which arbitrary elements of SU (2) are utilized. In such a set up, is it 
possible to use quantum game theory to come up with a natural choice for the stable 
state? Moreover, is it possible to characterize a quantized version of the Parrondo effect 
in this general set up, and if so, what does it mean for quantum computation? 
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To be more precise, the work in Chapter 3 embeds classical history dependent Par- 
rondo games into quantum multiplexers via embeddings of type 1 and 2. The resulting 
quantum multiplexers, when made to act upon a particular evaluative initial state, repro- 
duce the payoff functions of the classical Parrondo games. Call such quantum multi- 
plexers mundane. In other words, mundane quantum multiplexers reside in the image of 
the embeddings of type 1 or 2. However, the set of quantum multiplexers is much larger 
than the image of embeddings of either type; that is, there are quantum multiplexers that 
are outside such an image. Call such quantum multiplexers exotic. 

Clearly, the answer to question 1) above is in the affirmative for mundane quantum 
multiplexers based on the results of chapter 3. By taking quantum superpositions of the 
mundane quantum multiplexers associated with classical Parrondo games, the payoff 
function of the classical game can be reproduced by choosing a particular evaluative 
initial state such that the game is winning, even when the individual quantum games 
were losing with respect to appropriate evaluative initial states. For exotic quantum 
multiplexers, the answer is not clear cut since it is not known what an evaluative initial 
state for such a multiplexer should be. Therefore, a future study toward answering 
question 1) requires efforts into identifying such an appropriate initial state for exotic 
quantum multiplexers. In the context of quantum logic synthesis, how might an arbitrary 
quantum logic gate be synthesized via decomposition in a game theoretically meaningful 
way? Tha it, first assign a fixed number of qubits in the circuit to each player. Then, for 
an arbitrary quantum logic gate U, how might U be decomposed into sets of one qubit 
gates, one for each player, and an initial state choosen, such that a given game theoretic 
outcome might be realized? 
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Appendix A 



QUATERNIONS 

Complex numbers are extension of real numbers. This fact motivates us to view quater- 
nions as extension of the complex numbers, with the exception that the recipe for con- 
structing the conjugate of a complex number needs modification when one tries to follow 
it to construct the conjugate of a quaternion. This modification is such that the quater- 
nionic product is necessarily non-commutative and satisfies zj — jz for any complex 
number z and the quaternion j. 

A.l Complex Numbers 

The set of complex numbers is 

C = {oo + a\x : oq, oi G R and — — l} . 

Since complex numbers are just first degree polynomials, one defines binary operations 
of addition and multiplication on C via polynomial addition and multiplication respec- 
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lively. 

Addition : (oq + aix) + {bo + bix) = (ao + &o) + (cn + bi)x 
Multiplication : (ao + aix){bQ + bix) — aofeo + aibox + aobix + aibix^ 

The constraint x'^ — —1 provides multiplicative closure to C, yielding 

(ao + aia;)(6o + bix) = (ao^o - 01^*1) + (00^*1 + 01^*0)2^ 

The equation = — 1 has exactly two solutions, x = and x = —1/^. Setting 
X = i = leads to the conventional notation for the complex numbers 

C = {ao + ttii : Qq, ai G R and = — l} . 

The solutions i and —i are called imaginary numbers. This terminology gives rise to the 
notion of the real part ao and the imaginary part ai of the complex number ao + aii 
Note that since —i is also a solution to the equation x'^ — —1, there are complex numbers 
in C of the form 

Oo + Oi(— i) = ao — aii. 
The latter is called the conjugate of the complex number ao + aii, and one checks that 

(ao + aii)(ao — aii) = ag + a^ e M 

Clearly, the conjugate of ao — aii is the complex number ao + aii; that is, double 
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conjugation gives back the original complex number. The quantity 

\ao + aii\ = \jal + a\ 

defines the length of the complex number ao + oii (and of ao — oii). It is an easy exercise 
to show that C in fact forms a field. 

A.2 Quaternions 

The set of quaternions is 

H = {po : Po,Pi e C and = -l} . 

Again, addition and multiplication in EI is defined as polynomial addition and multipli- 
cation, giving 

Addition : (po + Piv) + (go + QiV) = {Po + Qo) + {Pi + qi)y 
Multiplication : {po + p^y) (go + Qiv) = Po^o + (Pi^o + Po^i)?/ + PiQiv'^ 

Is H closed under multiplication? The answer is yes once we note that Po, Pi? goi Qi are 
all complex numbers and that this requires the use of both the constraints = — 1 and 
— —I 'm simplifying the quatemionic product. Let 

Po = Po + P'lh Pi = Po + Pi^> % = % + Q'lh Qi = Qo + q'I^ 
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be complex numbers. Simplifying the quaternionic product now results in the expression 

(Po+Piy)(qo + qiy) 

= {pWo - pWi) + {pWi +PWo)i 

I r/ // / II I \ . I II I . II I \- . I I II I ll\ . / I II , I ll\ ■-{ 

+ [{PoQo-PiQi) + VPqQi+PiQq)^ + yPo% - PiQi) + [PoQi +Pi%)Ay 

+ [{pX-p'Wi) + {p'k'i+p'[W 

= {pWo - pWi - pWo - pWI) + (pWi + pWo - pWI - pWo)^ 

+ (pWo - pWi + pWo - pWDv + (pWi + pWo + PoQi + PWo)iy 

III II II II II ll \ , f I I , I I II II II ll\ ■ 

= [PoQo - PiQi - Po% - Pill) + {Po(li+Pi%-Po(li -PiQoP 

+ [{pWo - pWi + pWo - pWI) + m + pWo + pWI + pWo)i] y 

= Zo + ziy 
for complex numbers 

zo = {pWo - pWi - pWo - pWD + (pWi + pWo - PoQi ' pWo)^ 

and 

I II I II I , I II I ll\ , I II I , II I X ' " , I ll\ ■ 

zi = [Poqo-PiQi +Po% -PiQi) + [PoQi + PiQo + Poll +PiQoP- 

It is important to note here that even though the variable y is a square root of — 1, 
it is not equal to ±i. For if it were equal to ±i, then the set EI would equal the set C! 
By analogy with the complex numbers, the variable y might appropriately be called an 

imaginary complex number. It is commonly known as a hypercomplex number. FoUow- 
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ing convention, we replace y with j and write quaternions as po + Pij- 

We next develop the notion of a conjugate of a quaternion; that is, for a given quater- 
nion p, find a quaternion q such that pq E M. Following the recipe that led to the defini- 
tion of the complex conjugate naively we set po +Pi{—j) = Po~Pij as the quaternionic 
conjugate of the quaternion po + Pij- This gives 

(Po + Pij)(Po - Pij) =pI+ PiPoJ - PoPij +Pv (A.l) 

Multiplication of a complex number by its conjugate results in a real number that is 
the sum of the squares of two real numbers, namely the real and imaginary parts of the 
complex number. Since our definition of the quaternionic conjugate is motivated by the 
complex conjugate, we expect the right hand side of equation (A.l) to equal to the real 
number that results from the squares of the complex numbers po and pi . However, the 
fact that in general the square of a complex number is another complex number puts a 
kink in our plans. But all is not lost. Instead of insisting on the squares of the complex 
numbers po and pi in our definition of the quaternionic conjugate, we are perfectly 
happy to work with the squares of the lengths of the complex numbers po and pi, which 
are both real numbers. This flexibility forces us to modify the proposed quaternionic 
conjugate to the quaternion {p^ — pij) which gives 

(Po + Pij)(Po - Pij) = IpoI V p\p^3 - PoPij + bil^ (A.2) 
To eliminate the quaternionic part from the right hand side of equation (A.2) we are 



108 



Appendix A. Quaternions 



forced to set 

PiPo = PoPi = PiPo 

which means that pipo is in fact a real number, sacrificing the generality of our argument. 

At this stage, one wonders whether the recipe for the complex conjugate that has 
been followed thus far with a slight modification to develop the quatemionic conjugate 
needs to be changed drastically. Indeed, if we leave out the major ingredient of commu- 
tativity from the recipe and assume that for a complex number z, 

zj = jz, (A.3) 

then equation (A.l) must be re- written as 

(Po + Pij)iPo - Pij) =pI+ PijPo - PoPij + PijPij 
= pI+ PiPoj - PoPij + PiPijj 
= pI+ PiPoj - PoPiJ + bi 1^ 

The occurrence of \pif in the preceding equation is glaring, and suggests that we modify 
the proposed quatemionic conjugate yet again to be po — Pij which upon multiplication 
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with po + pij and after using the non-commutativity condition zj = jz leads to 



(Po + Pij)(po - Pij) = \pof + PijPo - PoPij + \p 



= \pof + piPoi - popij + bil^ 

= ip'oY + ip[r + ipif + iPif 

The quaternionic conjugate defined this way behaves much like the complex conjugate. 
For example, the quaternionic conjugate of po — pij is po + Pij- Moreover, as with 
the complex conjugate, the product of a quaternion with its conjugate is expressible as 
the sum of squares of four real numbers. We use the latter to define the lenght of a 
quaternion as 

bo + Pi3\ = Vbol' + bii' = vW+W+MF+MF- 

Rewriting po + pij as 

Po + Pij = {p'o + P'li) + {Po + KOi = Po + Pi^ + Poj + Pi^j (A-4) 

introduces the term ij which the non-commutativity condition of equation (A. 3) shows 
to be a square root of —1. For convenience, set k = ij. Then one computes 
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Complex number arithmetic together with equation (A. 3) establish the following identi- 
ties as well. 

ik = = i^j = -j 

jk = j{ij) = {-i)f = i 

The last two identities and the identity ij = k establish the right-hand rule for quater- 
nionic multiplication which is conveniently represented in the picture below. This rule 
is summed up in Hamilton's Relation P = = k'^ = ijk = —1. 

k j 




One can verify that the quaternions satisfy all the axioms of a field except commuta- 
tivity, and therefore form a division ring. Our definition of the quaternions in fact shows 
that the quaternions form a two dimensional algebra over the complex numbers with 
basis {1, j}. Equation (A.4) shows that the quaternions form a four dimensional algebra 
over the reals with basis {1, i, j, k}. 
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Appendix B 



COSINE SINE DECOMPOSITION OF UNITARY MATRICES 

As we shall see, the cosine sine decomposition (CSD) is essentially the well known 
singular value decomposition (S VD) of a unitary matrix implemented at the block matrix 
level. The reader is cautioned that for a given matrix, the CSD is not unique. The 
material presented in this appendix is not new. The discussion of the SVD is based on 
lecture notes of Professor Bin Jiang at Portland State University and the CSD discussion 
is based on the account given in [j-t] on pages 37-40. 

B.l Singular Value Decomposition 

Begin with the vector and matrix 2-norms, described below. 

Definition B.l. The 2-norm of a vector x G C" is the function || : C" M defined 



by 



X 



2 




Here, x"!" = (x^, X2, . . . , 



and |xjp = XjX* for Xj G C. 



Definition B.2. The 2-norm of a matrix A G C™^" is the function || \\^ : C 
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defined by 

\\Ax\\ 

\\-^\\2 

Since we will not refer to any other norms that can be defined on vectors and ma- 
trices, from now on we will use the || || instead of the more explicit || to simplify 
notation. Also, for A e C"^^", denote by A^ the conjugate transpose of A. Recall that 
a matrix A is unitary if AA^ — A^A — I. Equivalently, the action of a unitary matrix 
preserves vector norm. 

Lemma B.3. Vector and matrix 2-norms are invariant under unitary transformations. 
Proof. Let U e C"^^" be a unitary transformation, and x e C". Then 

\\Ux\\ = {{Ux)\Ux)Y = {x^U^UxY = {x'^xY = \\x\\ 

Now let A e C'"^". Then 

\\AU\\ = max||^||=i \\AUx\\ = max||[/^||=i \\AUx\\ = maxn^^n^i \\Ay\\ = \\A\\ 

If A e C"^". Then 

\\UA\\ = max||^||=i = maxn^^n^i = max||^||=i = \\A\\ 

□ 

We are now ready to prove the existence of a singular value decomposition. 
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Proposition B.4. If A e £rnxn^ ^^^^ there exists unitary matrices U e £rnxm 
V e C"^", and a matrix E = diag((7i, (72, . . . , dp, 0, . . . , 0) e R"^^", p = min(m, n), 
such that 

A = UEVl 

The (jj are called singular values of ^4 and are typically ordered so that 

Proof. The proof will be inductive. Let a — \\A\\. Since 

\\A\\ = max||^||=i \\Ax\\ , 

there exists a unit norm x e C" such that a — \\Ax\\; therefore. Ax — ay for some 
y e C"* with ||y|| = 1. 
If 

Vi ^{vi V2 ... Vr) e C"*^^ r<m 
has orthonormal columns Vi, then applying Gram-Schimdt process we can always find 

so that (V^i, V2) is unitary and rank(\4)"'" — rank(V2)- From this fact we conclude that 
there exist V{ e C"^^"-^) and ¥( e C'^xim-i) s^^h that Vi = {x ¥() e C"^" and 
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Ui = {y U[) e C"*^("^-^) are unitary. Hence, 



UlAVi 



A{x VI) 



^ y'^Ax y'^AVi ^ 
^{U'^^Ax {U'.YAVi ^ 

^ y^ay y^ AV[ ^ 
^{V[fay {UifAVl ^ 

B 



( 



\ 



= A^ 



) 



where w'^ e M^"-^). 

In fact w — For by lemma B.3., ||Ai|| = ||A|| = a and 



pill = maxy^ii^o 



||Aix| 



> 



\x\ 



/ \ 



J 

Bw 



> 



W 



w 
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Therefore, a > Vcr'^ + w'^w and hence w'^w — which implies that w — 0. 
We now have that 



UlAVi 



a 



(B.l) 



Now applying the same method to B and the resulting blocks B' inductively, we 

have 

C/J . . . C/|C/|^ 14^2 . . . l^p = diag((7i, (72, . . . (7p, 0, . . . , 0) 
Let \J — U1U2 . . .Up and VV1V2 . . .Vp. Then both U and V are unitary and 



□ 



B.2 Cosine Sine Decomposition 



Proposition B.5. Let the unitary matrix We C"^" be partitioned in 2 x 2 block form 
as 

I n — r 



W = 



n — I 



W21 W22 



I 



with 21 < n. Then there exist unitary matrices U — diag(t/ii, U22) and V = diag(Vii, V22) 
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with Un, Vn e C'^' such that 



I 
I 

n-2l 



I I n-2l 



c -s 
s c 








I 



(B.2) 



where 



C = diag(cos6'i, 0086*2, • • • , 0086*/) 



5" = diag(sin 9i, sin 6*2, ... , sin 6*/) 



such that sin^ 9i + cos^ 9i = 1 for some 9i, 1 < i < I. 



Proof. Let 



UlWiiVu = C 



be a singular value decomposition of the block H^n of and suppose that 



C = diag(Ci,/;_fe) 



where the diagonal elements of Ci satisfy 



< Ci < C2 < . . . Cfe < 1. 



Note that since W is unitary, the singular values cannot be greater than 1. Clearly, the 
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columns of the matrix 



11 



21 



are orthonormal. Therefore, 



Wii 
y W21 J 



n t 



y W21 j 



V^ii 



C + (W^2i^ii)^ (W^2i^ii) ; 



that is, 



(W^2i14i)^(W^2i1^ii) = diag(7 - Cl, 0^_fc) 



The columns of VF2i^ii are orthogonal with the last {I — k) of them being 0. Thus, there 
exists a unitary matrix U22 £ C^""'^^^""'^ such that 



where 



S = diag(si, S2, . . . , Sfc, 0, . . . , 0) = diag(S", 0) 



(B.3) 



with 5" consisting of k rows and the all O's block consisting of (r — A;) rows. Since 



diag(C/n,C/22)^ 



W^ii 



11 



w, 



21 



s 

\' J 
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has orthogonal columns, it follows that for 1 < i < I 

c^ + s^ = l. (B.4) 

In particular, S' is non-singular. 

Similarly, we may determine a unitary matrix V22 G £.i^-'-)>^(^~'-) such that 

UlW,2V22 = {T,0) 
where T = diag(ti, ^2, • • • , ti) with ti < 0. Since 

UliWn Wu)die.g{Vn,V22) = {C T 0) 

has orthogonal rows, it must be that cf + = 1, and it follows from (B.3) and (B.4)that 

T = -S. 

Now set U = diag(f/ii, f/22) and V = diag(Vii, V22)- Then it follows from the 
preceding discussion that 

X = U^WV 
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can be partitioned as 



X 







k 


/ - k 


k 


l-k 


n-2l 


k 









-Si 





\ 


l- 


k 





/ 











k 




Si 





-^33 




-^35 


l- 


k 










X44 




n - 


-21 


1° 







-^54 


-^55 y 



(B.5) 



Since X is unitary and Ei has positive diagonal elements, we have X33 = Ci. Moreover, 
X34, X35, X43, and X53 are zero. Therefore, the partition of X in (B.5) can now be 
updated to 







k 


l-k 


k 


l-k 


n-2l 


k 




' C 





-S' 





\ 


l- 


k 





I 











k 




S' 





c 








l- 


k 











X44 




n - 


- 21 










X^4 





(B.6) 



and the the matrix 



33 



( X X ^ 

A44 A45 

X54 X55 



is unitary. 
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Now we have 



diag(/('+'=),^l3)X 



Note that 



Hence, if we set 








-s. 











I 











^1 





Ci 

















I 

















/ 



I 



I I n-2l 



C -S 

s c 








n - 2/ y 



diag(7('+'=), Ul)X = diag(7('+'^), Ul)U'^WV. 



U = Udmg{I^'+^\Uss) 

= diag(C/ii, C/22)diag(7«, diag(7('=), C/33)) 
= diag([/ii,f/22-diag(7('=),t/33)) 
= diag(C/ii,C/22) 



Set 



C/2 = diag(7fe,C/3)C/2 
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and 

f/ = diag(t/i,t/2) 

Then 

[/tW = diag(/,+fc,f/3)X, 
then WWV has the form (4.2), where U and V are block diagonal unitary matrices. □ 
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LIST OF NOTATIONS AND NOMENCLATURE 

The state space of one qubit is the two dimensional complex projective Hilbert 
space CP^. As is the convention in quantum mechanics, an element i/j of the state 
space is denoted in Dirac notation by and is called a "ket" vector. 

|0) = (1,0)^ and |1) = (0, 1)^ are elements of the orthonormal computational 
basis of CP^. We point out that every ket is a column vector, however, as is the 
case here, it is sometimes written as the transpose of the appropriate row vector 
for notational convinience. 

1-0) = -00 |0) + -01 11) = ("00, V'l)^ is a quantum superposition of the elements of 
the computational basis, with iV'ol^ + iV'i 1^ — 1- In the language of linear algebra, 
lijj) is a linear combination of the elements of the computational basis. 

The dual of \ip) is the "bra" vector (^| = (^o '0i). where ijji is the complex 
conjugate of the complex number ^j. Note that a bra vector is a row vector. 

For l'^) = (■00, ipi)^ and |0) = (0o, 0i)^ in CP^, their inner product is given by 
(|'0) , 10)) = (■00 ■0i)(0o, 0i)^ and is denoted in the bra-ket notation by (■0| |0) 
or just (-010). 
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Appendix C. List of Notations and Nomenclature 

• The outer product of IV') and |0) is denoted by IV') (0| and is used to construct 
measurment operators. 

• If M is a matrix, then is the conjugate transpose of M. If M is unitary, then 

• The trace trace (A) of a square matrix A is the sum of its diagonal elements. 
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